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ABSTRACT 

Context. Type la Supemovae (SNe la) have been used as standardizable candles in the optical wavelengths to measure distances with 
an accuracy of ~ 7% out to redshift z ~ 1.5. There is evidence that in the near-infrared (NIR) wavelengths SNe la are even better 
standard candles, however, NIR observations are much more time-consuming. 

Aims. We aim to test whether the NIR peak magnitudes could be accurately estimated with only a single observation obtained close 
to maximum light, provided the time of B band maximum and the optical stretch parameter are known. 

Methods. We present multi-epoch UBVRI and single-epoch J and H photometric observations of 16 SNe la in the redshift range 
z - 0.037 - 0.183, doubling the leverage of the current SN la NIR Hubble diagram and the number of SNe beyond redshift 0.04. This 
sample was analyzed together with 102 NIR and 458 optical light curves (LCs) of normal SNe la from the literature. 

Results. The analysis of 45 NIR LCs with well-sampled first maximum shows that a single template accurately describes the LCs if 
its time axis is stretched with the optical stretch parameter. This allows us to estimate the peak NIR magnitudes of SNe with only few 
observations obtained within 10 days from S-band maximum. The NIR Hubble residuals show weak correlation with AM15 and the 
color excess E(B - U), and for the first time we report a dependence on the y^ax - Umax color. With these corrections the intrinsic 
NIR luminosity scatter of SNe la is estimated to be less than 0.08-0.10 mag, which is smaller than what can be derived for a similarly 
heterogeneous sample at optical wavelengths. We find that the intrinsic B - V color of SNe la at 6-band maximum light is a non¬ 
linear function of AM15 and estimate its intrinsic scatter to be ~ 0.04 mag. Analysis of both NIR and optical data show that the dust 
extinction in the host galaxies corresponds to a low Ry - 1.9. 

Conclusions. SNe la are at least as good standard candles in the NIR as in the optical and are potentially less affected by systematic 
uncertainties. We extended the NIR SN la Hubble diagram to its non-linear part at z ~ 0.2 and confirmed that it is feasible to 
accomplish this result with very modest sampling of the NIR LCs, if complemented by well-sampled optical LCs. With future facilities 
it will be possible to extend the NIR Hubble diagram beyond redshift z - I and our results suggest that the most efficient way to 
achieve this would be to obtain a single observation close to the NIR maximum. 
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1. Introduction 

Since the discovery of the accelerated expansion of the universe 
( |Riess et'ar]|1998| [Perlmutter et al.||l999| l the world sample 
of Type la supernovae (SNe la) at cosmological distances has 
grown substantially. Mainly thanks to the combined effort from 
the SDSS-II and SNLS surveys, the recent optical SN la Hubble 


Send offprint requests to\ vallery.stanishev@gmail.com 
* Partly based on observations made with ESO telescopes at the 
Paranal Observatory under program IDs 079.A-0I92 and 08 LA-0734. 


diagram was built from 740 well measured objects ( |Betoule et al. 
|2014| and references therein). These measurements are provid¬ 
ing the most accurate study of the expansion history of the uni¬ 
verse to date, especially when combined with CMB constraints. 
However, it is unclear at this point if the precision of the mea¬ 
surements of cosmological parameters, most notably the equa¬ 
tion of state parameter of the dark energy, will improve signifi¬ 
cantly from further increase in the sample of SN la at rest-frame 
optical wavelengths. Systematic uncertainties, instrumental, as- 
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trophysical and cosmological, may pose severe limitations (see 
Goobar & Leibundgut|201 1] for a recent review). 


One area of concern is the non-standard color-brightness re¬ 
lation observed in SNe la, both nearby and in cosmological sam¬ 
ples, implying that the current assumption of a universal redden¬ 
ing law independent of redshift may not hold. Similarly, the light 
curve (LC) shape-brightness relation used to calibrate SNe la 
lacks strong theoretical explanation, i.e., a redshift evolution 
cannot be excluded at this stage. 

One alternative that has been proposed is to explore the use 
of near-infrared (NIR; 1-2.5 pm) observations of SNe la. Ob¬ 
serving SN la in the NIR has several important advantages. First, 
the extinction corrections are by a factor of 4 - 6 smaller than 
in the optical B band. Second, there has been evidence that at 
maximum SNe la are natural standard candles in the NIR and no 


correction for the LC shape is needed for normal SNe (Meikle 

|2000| IKrisciunas et al.||2004a|cl Wood-Vasey et al. 

200811. Al- 

though more recent works by Folatelli et al.|(2010 

1 and |Kat-| 


tner et al. (2012 1 do suggest that such a correction may also 


be needed in the NIR, it is much smaller than in the optical. 
Third, |Krisciunas et al.| ( |2004c|b| l and later |Wood-Vasey et al.| 
( 2008| l demonstrated that smooth LC templates, the time-axis of 
which is stretched according to the SN optical stretch parame¬ 


ter Sopt (Perlmutter et al.|1997 1 , accurately describe rest-frame J 


and //-band LCs within 10 days from B-band maximum. This 
raises the possibility to measure accurately the peak NIR mag¬ 
nitudes of SNe la with even a single observation within a week 
from B-band maximum, provided the time of B-band maximum 
and iopt are known. This could be very important because due 
to the large atmospheric background, ground-based observations 
in the NIR rapidly become prohibitively time-consuming when 
observing faint objects and building well-sampled LCs of high- 
redshift SNe la becomes nearly impossible. The main goal of 
this work is to test the feasibility to measure the NIR peak mag¬ 
nitudes with single observation and extend the NIR Hubble dia¬ 
gram to intermediate redshifts. We report on the optical UBVRI 
and NIR J and H measurements of 16 SNe la in the redshift 
range z = 0.037 - 0.183. We present the results of the com¬ 
bined analysis of this new data-set together with a large set of 
published SN la light curves. 

The paper is organized as follows. In Sec. we present the 
target selection, observations, data reduction, photometry and 
calibration. Section [^provides details on the light curve htting, 
including the derivation of new NIR LC templates and NIR K- 
corrections. In Sec. we present our results and in Sec.j^a dis¬ 
cussion followed by our conclusions. Throughout the paper we 
assume the concordance cosmological model with = 0.27, 
Qa = 0.73 and h - 0.708. 


2. Observations and data reduction 

2.1. Target selection and observing strategy 


Between 2007 and 2008 we observed sixteen SNe la at redshifts 
between z = 0.037 and z = 0.183. Details of the targets are 
given in Table [1] The targets were selected from the SNe dis¬ 
covered by SDSS-II - 3 SNe, SNfactory ( Aldering et al.|20()2] l 
- 12 SNe, and ROTSE - 1 SN, to be spectroscopically classihed 
before maximum light. The early classihcation was essential in 
order to obtain the NIR observations within a week from max¬ 
imum and to obtain good enough optical observation to derive 
the LC parameters. For each SN, a single NIR observation and 
multiple UBVRI observations with cadence of 5 - 7 days were 
obtained, except for SNe 2007hz, 2007ih and 2007ik for which 
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the gri photometry obtained by SDSS-II was used ( |Sako et ah] 
2014| l. For most SNe the observations started earlier than 10 days 


after the B-band maximum light. 

Most of the SNe in our sample are close to their the host 
galaxy nuclei and the underlying non-uniform galaxy back¬ 
ground may signihcantly degrade the photometry. To remove the 
galaxy background contamination we used images of the host 
galaxies of the SNe that we obtained about a year after the main 
survey, when the SN had faded. In the NIR, template images 
were obtained for all three instruments used. In the optical, tem¬ 
plate images were only obtained with ALFOSC at NOT. The 
galaxy template images were aligned with the SN imag^ con¬ 
volved with a suitable kernel so that the point-spread functions 
(PSF) of the two images are the same, then scaled to match the 
flux level of the SN image and subtracted. The SN flux can then 
be correctly measured on the background-subtracted image. The 
image subtraction was done with Alard’s optimal image sub¬ 
traction software ( |Alard & Lupton|[l998[ Alard|[2000| , slightly 
modihed and kindly made available to us by B. Schmidt. SNe 
2007hz, SNF20080512-008 and SNF20080516-022 were found 
to be away of their host galaxies and the galaxy background in 
the NIR was not a problem. For these SNe no galaxy reference 
images were obtained and the photometry was performed with¬ 
out image subtraction. 


2.2. Near infrared imaging 

The NIR photometry of the SNe was obtained at three facili¬ 
ties. J and H imaging of the three SDSS-II SNe was obtained 
at VLT/ISAAC ( |Moorwood et al.||l998| l and nine of the others 
were observed with NOTCAM at the Nordic Optical Telescope. 
Six SNe were observed at the Japanese InfraRed Survey Facility 
(IRSF) 1.4m telescope located at SAAO and equipped with the 
imaging camera SIRIUS, which can simultaneously observe in 
JHKs bands. Two of the SNe were observed both at IRSF and 
NOT. All three instruments are equipped with 1024x1024 pixel 
Hawaii detectors and hlters from the Mauna Kea Observatories 
Near-Infrared Filter Set (MKO) ( |Tokunaga et al.|2002 1 . 

The observations consisted of series of dithered images in 
order to facilitate the night sky background subtraction. Twilight 
sky flat-held images were also obtained. The image reduction 
followed with small modifications the scheme outlined in IStan-l 
|ishev et aL| ( |2009| . All ISAAC images were first corrected for 
the effect of ’electrical ghost’ with the recipe ghost within the 
ESO software library Then dark current and flat-held 

corrections were applied to all images. The sky background sub¬ 
traction was performed with the XDIMSUM package in IRAFj^ 
To estimate the sky background for given image, the 8 exposures 
closest in time were scaled to the same mode and combined with 
median. The combined sky image was scaled to the mode of 
the image and subtracted from it. The pedestal on NIR arrays is 
highly variable in time and may be very difficult to remove ac¬ 
curately. In some of the images, significant residuals were found 
along the bottom and the middle rows of the array, where the 
readout of the two halves of the array starts. To correct for this 

' forth-order B-spline interpolation was used for all image interpola¬ 
tions in this paper because it has one of the best interpolation properties 
of all known interpolation schemes, e.g. IThevenaz et al.|P0001l. 

^ available at http: //www. eso. org/pro jects/aot/eciipse/ 

^ All data reduction and calibration was done in IRAF and with our 
own programs written in IDL. IRAF is distributed by the National Opti¬ 
cal Astronomy Observatories, which are operated by the Association of 
Universities for Research in Astronomy, Inc., under cooperative agree¬ 
ment with the National Science Foundation. 
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Table 1. Details of the SNe observed in this work. 


SN ID 

RA (J2000) 

Deck (J2000) 

■^helio 

ZCMB 

E(B - y)Mw 

Discovery 

SN2007hz 

21:03:08.95 

-01:01:45.1 

0.1393 

0.1381 

0.070 

SDSS-II 

SN2007ih 

21:33:10.77 

-00:57:36.5 

0.1710 

0.1697 

0.042 


SN2007ik 

22:38:53.72 

-01:10:01.5 

0.1830 

0.1815 

0.043 


SN2008bz 

12:38:57.74 

H-11:07:46.2 

0.0603 

0.0614 

0.025 

ROTSE 

SNF20080510-005 

13:38:24.29 

h-09:40:1L3 

0.0840 

0.0850 

0.024 

SNEactory 

SNF20080512-008 

17:03:35.14 

h-31:21:1L7 

0.0840 

0.0840 

0.030 


SNF20080512-010 

16:11:04.35 

h-52:27:09.9 

0.0630 

0.0631 

0.018 


SNF20080516-000 

11:22:29.71 

h-56:24:26.5 

0.0720 

0.0726 

0.009 


SNF20080516-022 

14:59:27.97 

h-58:00:09.8 

0.0740 

0.0742 

0.013 


SNF20080517-010 

12:36:23.88 

h-39:30:43.8 

0.1200 

0.1209 

0.014 


SNF20080522-000 

13:36:47.59 

h-05:08:30.4 

0.0450 

0.0460 

0.023 


SNF20080522-001 

12:51:39.88 

-11:59:23.7 

0.0490 

0.0502 

0.039 


SNF20080522-004 

13:11:32.70 

-18:39:45.3 

0.1050 

0.1062 

0.069 


SNF20080522-011 

15:19:58.91 

h-04:54:16.9 

0.0370 

0.0376 

0.037 


SNF20080606-012 

17:43:57.23 

h-51:22:18.8 

0.0790 

0.0788 

0.031 


SNF20080610-003 

16:17:37.73 

h-03:28:28.7 

0.0930 

0.0933 

0.054 



effect the median value of each row was computed with 2cr clip¬ 
ping and subtracted. The sky-subtracted images were combined 
and the combined image was used to create an object mask with 
the IRAF OBJMASK task. The mask was then de-registered to 
each individual image and the sky subtraction was repeated, but 
now all pixels belonging to objects were excluded during the 
background estimation and the residual pedestal subtraction. The 
only difference was that instead of the median, the mode along 
the image rows was used to remove the residual pedestal. 

Before the final image combination residual cosmic ray hits 
were identified and cleaned with the Laplacian detection algo¬ 
rithm of jvan Dokkum| ( [20OT| ). ISAAC and NOTCam have signif¬ 
icant optical image distortions. We derived corrections for this 
effect using observations of dense stellar fields with astrometry 
from SDSS and combined the images with the DRIZZLE algo¬ 
rithm ( [Fruchter & Hook1|2002[ ). For the final combination, the 
individual images were also assigned weights in order to opti¬ 
mize the signal-to-noise ratio (S/N) for faint point sources. The 
weights are inversely proportional to the product of the square 
of the seeing, the variance of the sky background and the multi¬ 
plicative factor that brings the images to the same flux scale. 


2.2.1. Photometry and calibration 

The instrumental magnitudes were measured by PSF photome¬ 
try with the DAOPHOT package ( |Stetson|ri987[ ) in IRAF. The 
photometric calibration of the SNe observed at ISAAC was ob¬ 
tained with standard stars from [Persson et al.| ( |I998| ), which 
were observed as part of ISAAC standard calibration plan. The 
zero-points estimated from the standards were transferred to 
the SN images assuming atmospheric extinction coefficients of 
kj = 0.05 and kn = 0.04 magairmass^*. The image zero points 
and their uncertainties for the observations obtained at NOT and 
SAAO were computed with 2MASS stars in the field-of-view. 
The final calibrated magnitudes are given in Table 


2.3. Optical imaging 

For SNe 2007hz, 2007ih and 2007ik the photometry obtained by 
SDSS-II was used (Sako et al. 2014| l. For the remaining 13 SNe, 
the optical observations were obtained with four different instru¬ 
ments, ALFOSC, MOSCA and STANCam at the Nordic Optical 


Telescope (NOT) equipped with broadband UBVRI filters and 
ANDICAM at CTIO 1.3m telescope operated by the SMARTS 
Consortium equipped with BVRI filters. Standard fields from 
|Landolt| ( |2009| l were also observed with each instrument. The 
observations were collected between May and July 2008. On 
May 23rd, 2009 the SN fields were re-observed with ALFOSC 
to obtain reference images for the galaxy subtraction. All CCD 
images were bias and flat-held corrected, and trimmed. Cosmic 
ray hits were identified and cleaned with the Laplacian detection 
algorithm of van Dokkum P001|l. When multiple exposures per 


Alter were obtained, they were combined into a single image. 


2.3.1. Photometry and calibration 


For the optical observations aperture photometry was chosen 
over PSF photometry because the main instrument used for the 
survey, ALFOSC, has significant PSF variation over the image 
and in most SN fields there were not enough bright stars to model 
this effect. The instrumental magnitudes were measured through 
digital apertures with diameter of 2x FWHM. 

The standard star fields observed in eight photometric nights 
with ALFOSC at NOT were used to calibrate local sequences of 
stars around the SNe. Aperture corrections from 2xFWHM to 
5 X FWHM were computed with IRAF.PHOTCAL package and 
applied to the measured magnitudes. All standard star measure¬ 
ments were fitted simultaneously (with 3cr clipping) with linear 
equations of the form ( [Harris et al.|198l| : 

U = u + ctu{U-B)-k’jjX-k'{jX(U-B)+zpu 
B = b + ctB(B-V)-k'^X-k”X{B-V) + zpB 
V - V + ctv (B - V) - kyX + zpv (1) 

R - r + ct/itV - R) - k^X + zpR 
I — i + ctj (V - I) - k'j X + zpi, 


where upper-case and lower-case letters denote the standard and 
instrumental magnitudes, respectively, ct the color-terms, X the 
airmass, k' and k" the first- and the second-order extinction coef¬ 
ficients and the zp the instrument zero points. The second-order 
extinction coefficients k" were fixed to -0.03 based on our syn¬ 
thetic photometry calculations and past photometric works (see 
also Bond|20()5] l. Because the observations were obtained during 
a short period of time the data were fitted with common color- 
terms and instrument zero-points. The standard star observations 
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Table 2. NIR supernova photometry. 


SN ID 

JD-2400000 

J 

H 

Phase" 

Instrument 

sn2007hz 

54359.59 

20.904 (0.056) 

20.974 (0.096) 

7.2 

ISAAC 

sn2007ih 

54361.63 

21.641 (0.086) 

- 

6.5 


sn2007ih 

54359.67 

- 

21.587 (0.184) 

5.1 


sn2007ik 

54362.55 

- 

21.444(0.226) 

4.5 


SNE20080510-005 

54609.48 

20.215 (0.195) 

19.919(0.171) 

4.0 

NOTCAM 

SNE20080512-008 

54609.55 

20.579 (0.095) 

20.037 (0.102) 

13.7 


SNE20080512-010 

54608.55 

19.459 (0.044) 

19.475 (0.064) 

7.2 


SNE20080516-000 

54608.50 

19.077 (0.051) 

19.222 (0.067) 

1.9 


SNE20080516-022 

54609.50 

19.200 (0.045) 

19.213 (0.084) 

-0.3 


SNE20080517-010 

54610.46 

20.697 (0.102) 

20.888 (0.270) 

4.0 


SNE20080522-000 

54627.46 

18.319 (0.033) 

18.190 (0.031) 

5.5 


SNE20080606-012 

54633.47 

19.183 (0.027) 

19.461 (0.086) 

1.5 


SNE20080610-003 

54633.54 

19.890 (0.051) 

19.770 (0.080) 

0.1 


sn2008bz 

54592.31 

20.023 (0.265) 

18.780 (0.210) 

13.0 

SIRIUS 

SNE20080510-005 

54609.30 

20.063 (0.141) 

19.564 (0.150) 

3.8 


SNE20080522-000 

54614.28 

17.703 (0.018) 

17.930 (0.036) 

-6.0 


SNE20080522-001 

54613.31 

18.965 (0.062) 

18.719(0.091) 

8.8 


SNE20080522-004 

54614.46 

19.861 (0.130) 

19.988 (0.276) 

5.8 


SNE20080522-011 

54614.37 

17.284 (0.017) 

17.547 (0.036) 

-1.9 



Notes. stretch-corrected rest-frame phase (Eq|^. 


with the other three instruments we only used to derive the color- 
terms, again assuming k" - -0.03 for U and B bands. 

With the derived transformation equations, the standard 
magnitudes of the stars in the SN fields were computed. The 
weighted mean and standard deviation of the calibrations from 
different nights (weighted by the estimated photometric uncer¬ 
tainties) were taken as standard magnitude and its uncertainty. 
The stars that showed scatter larger than that expected from the 
photometric errors were discarded. The final calibrated local se¬ 
quences are indicated on the finding charts in Fig.[T]and the mag¬ 
nitudes are listed in Tabled 

The calibrated local sequence stars, measured with 
2xFWHM diameter apertures, were used to compute the image 
zero points. For each SN image, a zero-point zp^ = zp-k' X was 
calculated for each calibrated star by applying Eqs.[T] The final 
image ZPs and their uncertainty are, respectively, the average of 
the individual ZPs (with 3cr outliers removed if present) and the 
standard deviation. The image ZPs were added to the measured 
SN magnitudes measured in the galaxy-subtracted images to ob¬ 
tain the SN magnitudes in the natural systems of the instruments 
used, m™*. 

It is known that Eqs.[2cannot accurately transform SN pho¬ 
tometry into the standard Johnson-Cousins UBVRI photometric 
system ( Suntzeff||2000[ [Stritzinger et al.|2002[ Krisciunas et al. 
|2003| l. Eor this reason the optical photometry was calibrated with 
the method described in |Stanishev et al. ( 2007|l, wh ich is based 
on the 5-correction introduced by Stritzinger et al. (2002i. Be¬ 
cause we take into account the second-order extinction, the cal¬ 
ibration equations were slightly modified. The calibration cor¬ 
rection from natural magnitudes at airmass X, in'^\ to standard 
magnitudes is computed by means of synthetic photometry: 


m = m 


-2.5 log (/ fJ^°\A)R^^\A)dA) 

-H2.51og(/ ff°\A)R''’^\A)p(A)^dA) +const. (2) 

Here /]’*’°'(/l) is the photon spectral energy distributions (SED) 
of the SN, f?'’‘‘'(/l) and 7?®'‘*(/l) are the responses of the natural 
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and the standard systems, respectively. Note that f?™‘(/l) does not 
include the atmospheric transmission and thus the transmission 
of the atmosphere at airmass X - 1 above the site of observation, 
piA), is also included in the calculation. Because the SNe in our 
sample do n ot have extensive spectrophotometry we used|Hsiao| 


et al. 


(|2007| spectral template as /J^°’^(/l), corrected for the Milky 
Way reddening and red-shifted to the redshift of each SN. The 
constant in Eq.|^is such that the correction is zero for AO V stars 
with all color indices zero. This ensures that for normal stars 
the synthetic 5-correction gives the same results as the linear 
color-term corrections (Eq. [T]). The constant can be determined 
from synthetic photometry of stars for which both photometry 
and spectrophotometry is available. In the Appendix we provide 
full account of the calibration procedure, including a discussion 
of the uncertainties introduced when the second-order extinction 
is not taken into account. The final photometry of the SNe is 
given in Table 


3. Light curve fitting 

3.1. Template light curve fitting 

There are several different EC fitting methods available, stretch 



that works in both optical and NIR. However, we chose to de¬ 
velop and use our own simpler stretch-based fitter for the follow¬ 
ing reasons. Eirst, our primary goal is to test the accuracy of the 
estimate of SN la NIR peak magnitudes from a single observa¬ 
tion obtained close to maximum. |Krisciunas et al.] ( |2004b] l have 
presented evidence that the stretch approach could be applied to 
the first maximum of J band LCs and that the NIR and the opti- 
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sn2007hz z=0.139 s=1.069 sn2007ih z=0.171 s=1.184 sn2007ik z=0.183 s=1.090 


sn2008bz z=0.060 s=0.935 



0 10 20 30 40 50 

SNF20080522-004 z=0.105 s=1.249 


0 10 20 30 40 

SNF20080522-011 z=0.037 s=1.052 


-10 0 10 20 30 40 50 

SNF20080606-012 z=0.079 s=0.954 


10 20 30 40 50 60 70 

SNF20080610-003 z=0.093 s=1.047 



t-k 




Wr. 


i-Ir. 


Fig. 2. griJH and UBVRIJH LCs and the best fits. For the three SDSS-II SNe (2007hz, 2007ih, 2007ik) from bottom to top are plotted g, r- 0.4, 
- 1, 7 - 2.8 and 77 - 3.8 (no J for SN 2007ik). For the other 13 SNe are plotted U + 0.8, B,V- 0.2, R - 0.6, 7 - 1.6, 7 - 3 and 77 - 3.5. Note that 
the simultaneous gri or UBVRI fits were only used for these plots. In the analysis we use the results from the simultaneous fits of the g and r, and 
B and V bands only. 


cal stretch are tightly correlated. With the stretch fitting approach 
we can now verify this result with a much larger sample of SNe 
with well-sampled LCs. Second, we use our own prescription to 
compute NIR ^T-corrections. It is based on a collection of NIR 
spectrophotometry of nearby SNe la, which also have NIR pho¬ 
tometry. Each spectrum was carefully checked and when needed 
corrected to match the photometry. SNooPy ( [Bums et al.|201l| 
uses the NIR part of|Hsiao et al.|(|2007|) spectral template, which 


was warped to match the mean color evolution inferred from sev¬ 
eral SNe la. We find that this template gives slightly different K- 
corrections from ours. Third, from the optical LCs we need 
B and V peak magnitudes, and a LC shape parameter. Consid¬ 
ering the above, stretch-based fitter appears to be best suited for 
our purposes. 
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Table 5. Stretch-AMi 5 correspondence. 


Stretch 

0.65 

0.67 

0.70 

0.75 

0.80 

0.85 

0.90 

0.95 

1.00 

1.05 

1.10 

1.15 

1.20 

1.34 

1.47 

AMi5 

2.000 

1.926 

1.835 

1.653 

1.512 

1.353 

1.212 

1.145 

1.034 

0.937 

0.893 

0.813 

0.766 

0.690 

0.516 


Our fitting program is based on the IDL fitting package MP- 
fiiQ It can fit multi-color LCs simultaneously using two ap¬ 
proaches. The first is to fit the data with rest-frame template LCs 
coupled with a spectral template to compute ^T-corrections ( |Kim 
|et al.|[T996] l in order to transform from rest-frame to observed 
magnitudes. In the second approach a spectral template is red- 
shifted and directly integrated through the observed filters. The 
latter approach requires that the spectral template is calibrated 
to match the rest-frame LCs of SNe la, which is the case for 
|Hsiao et^ ( |2007] l template. This template is generated in such 
a way that it matches the UBVRI light curve templates of Knop 
|et al.| ( |200^ and its B - V color at maximum is -0.058 mag. In 
both approaches, the fitted parameters are the common stretch 
s, the time of B-band maximum light and the magnitude 
offsets needed to fit each SN light curve. The LC templates are 
normalized to zero magnitude at and because each template 
magnitude offset is fitted independently, the procedure returns 
the rest-frame SN magnitudes at The program is flexible 
and one or several fit parameters can be fixed. It also allows the 
spectral template to be reddened by arbitrary reddening laws. 

The stretch parametrization can only be used for U, B or V 
bands. In the redder bands SNe la show a second maximum (or 
just a shoulder in R) and the stretch approach can only be applied 
to the first maximum (e.g., |Krisciunas et al.|2004b| . Therefore, 
the fits were limited to the observations obtained up to phase 
-1-10 stretch-corrected days in J and H bands, and up to phase 
-t-35 days in the optical. The stretch-corrected SN phase, (p, is 
defined as: 


= C 




i(l H-z)’ 


(3) 


^ the template was modified using Fitzpatrick (1999 i extinction law, 
allowing for negative reddening when needed. 


redshift of the filters effective wavelength caused by reddening. 
In particular, the derived stretch parameters are largely free from 
this effect. The optical LC parameters of the SNe from our new 
sample are shown in Table They indicate that there are no un- 
derluminous or highly-reddened SNe in our sample. Further, the 
optical observations of all but two SNe started before phase -1-15 
days and the NIR observations of only two SNe were obtained 
after phase -i-10 days. For illustration purposes, the simultaneous 
gri or UBVRI fits are shown in Fig.|^ In the analysis we use the 
results from the simultaneous fits of the g and r, and B and V 
bands only. 

The optical LCs for the literature NIR sample were taken 
from the following sources: Calan-Tololo (|Hamuy et al.|19961l, 
CfAl , 2,3 an d 4 ( |Riess et al.|1999[|Jha et al.|2006t|Hicken et al.| 
2009 |2012| l, Berkeley ( |Ganeshalingam et al. 2010| l, Carnegie 


et al.pOOSll, K. Krisciunas and collaborators ([Krisciunas et al. 


Supernova Project (CSP) ([Contreras et al.||2010[ [Stritzinger 


et al.||2011|l. Supernova Cosmology Project (SCP) ([Kowalski 


2004c|b||2603|l2001l|2000|[2007ll, SN 199 8bu ( [Hernandez etak 


2000l|Jha et al.|1999l|Suntzeff et al.|1999|l, SN 2000E (jValentini 
et al.ib03|l, SN 2002W ^Benetti et al.|2004ll, SN 20()2dj (jPig 


(Stanishev et al. 


torello et al. 2007 


nata et al. 2008 1 , SN 2003cg (jElias-Rosa et al. 2006|l, SN 2003du 
pOOTl l, SN 2005cf ( jWang et al.|[2009bl jP^ 


tj), SN 2008Q (Stanishev et al., in preparation). 


SN 201 Ife ( Munari et al. 2013]l. E or the SNe observed in the 


NIR by Barone-Nugent et al. (2012 1 (hereafter referred as 


BN12 


sample) the optical g and r LCs from Maguire et al. (2012 i were 


used. The CSP photometry is given in the natural photometric 
systems of the CSP telescopes and the observed B and V bands 
used were those fromjStritzinger et al.|(|201 1 1 . Eor the other SNe 


the standard Bessell B and V bands (jBesse 


1 & Murphy|[20T2ll 


where z is the SN host galaxy redshift, s is the stretch and the 
constant C is set to 1, unless specified otherwise. The stretch pa¬ 
rameter s can be either the optical stretch iopt or the NIR stretch 
snir, which will be defined later. 

3.2. Optical light curve parameters 

Eor the analysis of the NIR LCs we need the optical stretch 
■Sopt, B ^ peak magnitudes. To determine them, 

the Milky Way extinction-corrected g and r LCs of SNe 2007hz, 
2007ih and 2007ik by SDSS-II, and the B and V-band LCs of 
the remaining objects up to phase -1-35 days were fitted simul- 
taneousl}p] in a two-iteration procedure. In the first iteration the 


Hsiao et al. ( 2007 1 ) spectral template was integrated through the 
relevant filters and fitted to the data to obtain the SN B - V color 
^Smax’ “ ^)Bmax' ^hc template was then modified to match 
the observed SN (B - V)b„^ coloi^and the fitting was repeated 
to obtain the final parameters. This two-step procedure has the 
advantage that it accounts for the LC shape modification by the 


were used. If a given SN had optical LCs available from differ¬ 
ent sources, the optical parameters associated to it were derived 
from the LCs that had better phase coverage and sampling. In a 
few cases it was necessary to combine the optical photometry in 
order to obtain more reliable LC fits. 

Although analysis of the optical LCs is not the main goal 
of this paper, all SNe with LCs that allow stretch fitting were 
also analyzed. This analysis was performed because it may be 
relevant and of practical interest for SN cosmology. Eor many 
current and future high-redshift SN samples only the rest-frame 
B and V band are in the optical window. The redder bands are 
red-shifted in the NIR and are/will be unavailable. The UV part 
of the spectrum is red-shifted to the optical, but the UV prop¬ 
erties of SNe la are much less understood and the peak magni¬ 
tudes appear to have larger intrinsic scatter ( [Brown et al.|2010] 
Milne et al. 2010 2013[). Given this, it is interesting to study 


available at http://cow.physics.wisc.edu/~craigm/idl/ 
idl.html 

We would like to stress here that in our definition iop, is the common 
stretch that makes the template fit both B and F-band LCs simultane¬ 
ously _ 


how well the SN la standard candle can be calibrated when only 
rest-frame B and V observations are available, now with much 
larger sample of nearby SNe. In addition, with this large sample 
of high-quality LCs we can study the SN intrinsic B - V color 
at tB^^, (B - V)b„„„, which is important for estimating the dust 
extinction in the SN host galaxies. 

Using a sub-sample of ~ 200 SNe with well-sampled LCs 
and observations within two days from tB„„ we derived a new 
iopt - AMi 5 relation for normal SNe la, where AM 15 is the B- 
band magnitude decline during the first 15 days after the time of 
B-band maximum light. The observed B-band LC were corrected 
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Table 6. Average difference between the B and V peak magnitudes, and 
AMi 5 for the SNe in common between the large data sets. Only the 
sets with more than 10 SNe in common are shown. The number in the 
parentheses is the Icr scatter units of 0.001 mag. 


Difference 

Nsn 

B 


y 

AMis 

Berkeley-CSP 

16 

0.031 ( 

:38) 

0.045 (25) 

-0.006 (44) 

Berkeley-CfA2 

10 

-0.027 

(38) 

-0.023 (14) 

0.032 (44) 

Berkeley-CfA3 

47 

-0.008 

(43) 

0.008 (45) 

-0.014 (54) 

CSP-CfA3 

30 

-0.050 

(35) 

-0.050 (24) 

-0.007 (56) 


to rest-frame B and normalized to zero peak magnitude. The LCs 
were divided into 0.05 stretch bins. Let be the stretch in the 
middle of a given bin. Then the time axes of the LCs that fall 
into this bin were converted to phases, corrected to the stretch in 
the middle of the bin by using Eq. with C - and stacked 
together. Smoothing spline was then used to estimate the time 
maximum and AM 15 parameter. The resulting Sopt - AM 15 pairs 
are given in Table From Fig. it is seen that our relation is 
close to the one derived by simply stretching the Z?-band LC tem¬ 
plate, indicating that the stretched template describes the LC of 
normal SNe la well. On the other hand the linear relation derived 
by Jha et al. ( |2006 1 , which has been used by many authors, is ap¬ 
plicable only in the stretch range 0.85-1.05. In our analysis we 
found that the "luminosity - light curve shape" relation for nor¬ 
mal SNe la is linear (within the scatter) when AM 15 is used as 
a LC shape parameter and not the stretch. Thus, we used AM 15 , 
which was computed from the stretch by cubic spline interpola¬ 
tion of the data in Table |5] 


Table 7. NIR spectra of SNe la used for computing Ai-corrections. The 
phases of the spectra are relative to the time of the B band maximum 
light. The sources of the NIR photometry are given in the text. 


SN 

Phase 

“^opt 

Ref. 

SN 1998bu 

-3,15,22,36 

0.973 

1,2 

SN 1999ee 

-9,0,2,5,8,15,20,23,32,42 

1.084 

2 

SN 2000ca 

2.7,9.7,17.7,22.6,30.8,38.6 

1.065 

3 

SN 2002bo 

11 

0.928 

4 

SN 2002dj 

-5.9 

0.948 

5 

SN 2003cg 

-6.5,-4.1,-3.5,4,14,19.6 

24,30.5,40 

0.937 

6 

SN 2003du 

-ll,-lL5,-5.5,-2.5,2.3,3.4 

4.5,12.4,16.2,20.4,30.4 

1.018 

7 

SN 2004S 

15 

0.992 

8 

SN 2005cf 

-11,-9,-6,-1,1,10,32,42 

0.979 

9 

SN 201 Ife 

-16,-14,-ll,-8,-4,-l,2,7,ll,16 

0.905 

10 


References. (l)|Jha et al.H1999J; (2) 

Hamuy et al.|(|2002j; (3) Stani- 

shev et al., in prep.; (4)|Benetti et al. i 

2004|l; (5) 

Pignata et al.|P008 ; 

(6) |Elias-Rosa et al. 1^2006 

; (7)|Stanishev et al.| 

2007|l; (8)|Gall et al.| 

(|2012|(; (9)|Pastorello et al. 

12007bll; (10)|Hsiao et al.H2013J. 


3.3. NIR K-corrections 

To compute the NIR Af-corrections we used spectrophotometry 
of nearby SNe la that also had NIR photometry available (Ta¬ 
ble]^. For the majority of the available spectra the J,H and K 
parts were obtained separately with different, non-overlapping 
instrument settings. While the relative flux calibration of each 
individual part is probably accurate enough, the cross-band flux 
calibration may not be. Besides, the flux across the strong telluric 
absorption bands at ~ 1.38 yum and ~ 1.88 /rm was not available. 
Therefore, before computing the AT-corrections the spectra had 
to be "prepared". The following procedure was applied. 

The noisy parts at the edges of the telluric bands were re¬ 
moved and the flux across the bands was linearly interpolated. 
Next, JHK synthetic magnitudes were computed with the ap¬ 
propriate effective system response^ and were compared with 
the observed ones. When necessary the J, H and K parts were 
scaled, the interpolation across the telluric bands was repeated 
and new synthetic photometry was computecj^ The procedure 
was iterated until the synthetic and the observed colors matched. 
In the cases when the whole NIR spectral range was observed 
simultaneously we always found a good match between the syn¬ 
thetic and the observed colors, and only interpolation across the 
telluric bands was necessary. The spectra were corrected for 
dust reddening in the Milky Way and the SN host galaxies us¬ 
ing E{B - V) values from the Galactic dust map ( |Schlegel et al. 
|1998| l and from the corresponding photometric papers, and cor¬ 
rected to rest-frame wavelengths. The spectra at similar phases 
were plotted together and examined. Several spectra with poor 
S/N, inadequate interpolation across the telluric bands or other 
problems most likely related to the flux calibration were identi- 


In Table we report the weighted mean difference and Icr 
scatter between the B and V peak magnitudes, and AM 15 of the 
SNe in common between the large data sets. Only the sets with 
more than 10 SNe in common are considered. One can see that 
there are small systematic differences in the peak magnitudes, 
which are significant to about Icr level. The reason for these 
systematic differences are likely related to the photometric cali¬ 
bration and different LC sampling. AM 15 shows relatively large 
scatter, but no significant bias. 


^ The effective system responses were computed by multiplying the fil¬ 
ter transmission, the detector quantum efficiency (QE), two aluminum 
reflections and the atmosphere transmission at airmass one. When avail¬ 
able the filter transmissions and the detector QEs were obtained form 
the instrument web-pages, if not, generic filter transmissions and QEs 
appropriate for the used instrument were used. 

® this is needed because in some cases the filters extend over the in¬ 
terpolated parts of the spectrum. As the J, H and K parts were scaled, 
the interpolation changes, which can also affect the computed synthetic 
magnitudes to some extent. 
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Fig. 4. Example of computing the NIR J and H Ai-corrections at z = 
0.02 (left panels) and z = 0.07 (right panels). 


fied and discarded. Details of the final set of spectra are given in 
Table [7] Note that a few spectra do not cover all JHK bands. In 
addition, few spectra had issues in one band, but were of suffi¬ 
cient quality in the others. 


To compute AT-corrections at any SN phase smoothing 
splines were used to interpolate the AT-corrections computed 
from the individual spectra. The smoothing splines have great 
flexibility to describe the data to any pre-defined rms around the 
ht. We have computed J- J and H-H AT-corrections for several 
redshifts between z = 0.002 and z — 0.2 and visually exam¬ 
ined the interpolation to tabulate the optimal smoothing param¬ 
eter as a function of the redshift. This is a somehow subjective 
procedure and we tried to select the smoothing parameters in 
such a way that the data was well-htted, but with special care 
to avoid over-fitting. Examples of the outcome of this procedure 
are shown in Fig. for two redshifts of z = 0.02 and z = 0.07. 
The scatter around the smooth curves increases with the redshift 
and also depends on partucular filters use (see below). The scat¬ 
ter at redshift z = 0.02 is ~ 0.01 mag in all bands, at z = 0.08 
it increases to ~ 0.033 and 0.04 mag, for J and H, respectively, 
and at z = 0.14 the J and H scatter is ~ 0.05 and 0.08 mag. 
These values and the observed increase with the redshift are sim¬ 
ilar to the results of |Boldt et al.| ( |20T4l l. Note, however, that we 
did not color-correct the spectra to follow singe color curve as 
in |Hsiao et al.| ( |2007| l and |Boldt et'af] ( |2014[ ). As a compari- 
son in Fig. |4 are also shown the AT-corrections computed from 


Hsiao et 2007 1 template. One can see that there are consider¬ 


able differences for some SN phases, including before -I-IO days, 
which is interest for our study. jBoldt et al.| ( |2014[ ) argue that the 
dominating part of the scatter is due to intrinsic differences in 
the energy distribution of different SNe. Thus, the scatter of the 
individual Af-corrections around the smoothing spline function 
between phases -10 and -I-IO days was tabulated as a function 
redshift and added in quadrature to the uncertainties of the peak 
magnitudes. 


There is no standard photometric system in the NIR and of¬ 
ten slightly different hlter sets are employed by different obser¬ 
vatories. In our analysis we decided to apply Af-corrections from 
the different observed filter systems to a single JHKs photomet¬ 
ric system, which we chose to be the one dehned by the Mauna 
Kea Observatories Near-Infrared Filter Set (MKO) ( |Tokunaga 
|et al.||2002] l. As a practical realization of the MKO system we 
selected the hlter/detector combination used by the NOTCam in¬ 
strument at the NOT. The observed filter systems that were used 
in this work are; 
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(i) the MKO filters for the 16 new SNe from this work and the 


SNe from Barone-Nugent et al. (20121 (BN12 set). 

(ii) 2MASS dCohen et al.||200^ or the CfA data set ( |Wood- 
Vasey et al.|20'08 1 


(iii) the filters used by The Carnegie Supernova Project (CSP), 
which are similar to |Persson et al.| (|1998|l. The CS P fil¬ 
ters were used for the CSP data set ( [Contreras et al.|[2010t 
Stritzinger et al.|2011 1 , the data set from K. Krisciunas and 


collaborators ( jKrisciunas et al.||2004c|b| [2003 1 [200 Ij |2000| 

|2007| ) (the Kri set) and the following SNe from various 
sources (the Var set) SN 1998bu ( [Hernandez et al. 2000[ 


Mayya et al.|1998|[Jha et al.|1999[l, SN 1999ee ([Hamuy et al. 
20021 SN 2000E ^Valentini et al.|2003[l, SN 2ote ^~ pgnata 


et al.|2008[), SN 2603cg ( [Elias-Rosa et al.|2006[), SN 2003du 
( [Stanishev et al.[[2007 1, SN 2003hv ( [Feloudas et al.[[2009[l. 


SN 2004eo ([Pastorello et al.|2007a| l, SN 2005cf ( [Wang et ai.[ 
[2009b[|Pastorello et al.|2007b[ l, and SN 2008Q (Stanishev et 
al., in preparation). 

(iv) the filters of WHIRC instrument were used for SN 201 Ife 
( [Matheson et al.|2012[ ). 

Strictly speaking the above procedure will work accurately 
only if the photometry is given in the natural photometric sys¬ 
tem defined by the filters. This is the case only for the CSP and 
CfA data sets, and to a large extent for the Kri data set. The ob¬ 
servations of the other SNe were obtained with filters that differ 
from those of the photometric systems to which the photometry 
was calibrated. Most of our new photometry is tied to 2MASS 
system. The rest of the SN photometry is tied to [Persson et al.[ 


([1998 1 or Hunt et al. (19981 systems, which are similar to each 


other. Because of this mismatch, it is expected that the photom¬ 
etry will contain small systematic uncertainties. To estimate the 
magnitude of this uncertainty we computed synthetic photome¬ 
try with the filters of the four photometric systems considered 
in this work at two epochs from the Hsiao spectral template - 
at NIR maxima and a week after We found that the er¬ 
rors should be smaller than 0.02 mag. The only exception was 
H band a week after In this case the error could be as big 
as 0.06 mag for most SNe and only for WHIRC observations of 
SN 201 Ife it could reach 0.1 mag. 


3.4. NIR template light curves 

There are currently several NIR FC templates available but with 
excellent new NIR photometry of many SNe la published re¬ 
cently (e.g., [Contreras et al.[|2OT0t [Stritzinger et al.[|2011[ ), we 
decided to derive new NIR LC templates. From the sample of 
SNe with both NIR and optical photometry we discarded known 
peculiar objects (SN 2002cx-likes, SNe 2006bt, 2006ot, and SN 
1991bg-likes). Further, only the SNe that are suitable for deriv¬ 
ing NIR FC templates were selected, e.g. SNe with several ob¬ 
servations obtained before phase -I-IO days and well-distributed 
along the FC. 

We started by fitting the J and //-band FCs simultaneously 
using [Wood- Vasey et al.| ( |2008[ l templates and our new NIR K- 
correction procedure. The templates were stretched with the op¬ 
tical stretch iopt and the zero time was fixed to as in|Krisci- 
[unas et aL] ( [20()4c|b| ). Examination of the well-sampled LCs re- 
vealed that some fits were not very good. The reason appeared to 
be that the optical stretch Sopt did not provide adequate stretching 
to fit the NIR FCs. To further test this we selected a sub-sample 
of SNe with well-sampled FCs, which we call "GOFD" sample. 
For a SN to enter this sample it has to have (i) at least 5 points in 
J and/or H band obtained before phase -I-IO days; (ii) the points 
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Fig. 5. The data used to create the J and H light curve templates. The 
templates are overplotted with red lines. (T„ is the weighted RMS of the 
data around the template during the phases earlier than -l-lO days. 


are evenly-distributed along the LC; (iii) has observations within 
2 days from the time of NIR maximum. Such SNe allow to per¬ 
form the fit of the first NIR LC maximum with the NIR stretch 
iNiR and the time of maximum as free fitting parameters. The 
40 objects in the GOLD sample are: SN 1998bu, SN 1999ee, 
SN1999ek, SN2000E, SN2001cz, SN2001el, SN2002bo, 
SN2003du, SN2004eo, SN2004ey, SN2005M, 
SN2005el, SN2005eq, SN2005hc, SN2005iq, 
SN2005ki, SN2006D, SN2006X, SN2006ax, 
SN2006kf, SN20061e, SN20061f, SN2006os, 
SN2007S, SN2007af, SN2007bc, SN2007bm, 


SN2002dj, 

SN2005cf, 

SN2005kc, 

SN2006et, 

SN2007A, 

SN2007ca, 


SN20071e, SN2008Q, SN2008bc, SN2008fp, 
SN2008gp, SN2008hv and SN2011fe. SNe 2005el, 2005eq, 
2005iq, 2006X and 2006ax have LCs from both CSP and CfA, 
and the LCs were analyzed separately. The fit of the CfA LC 
of SN 2005iq and SN 2006ac provided very weak constrains on 
the stretch and the time of maximum and was moved from the 
GOLD sample to the SILVER one (see Sec. 3.5 1 . The differences 
between the peak magnitudes estimated from the CSP and CfA 
photometry are at most 0.08 mag. The only exception is H band 
of 2006ax, for with the //-band photometry from CfA was 0.22 
mag fainter than CSP. Even though the analysis of the maximum 
colors indicated that those obtained from the CSP photometry 
were more likely to be correct also kept the CfA photemetry of 
this SN. 


To derive the templates the observed LCs were X-corrected 
to rest frame magnitudes, normalized to zero magnitude at 
converted to phases using Eq. [^with s - snir, and stacked to¬ 
gether. The data for phases before H-10 days were fitted with 
smoothing-spline function, manually adjusting the smoothing 
parameter to achieve a smooth fit that leaves no correlated resid¬ 
uals and again taking care to avoid overfitting. A few points that 



^opt 



Fig. 6. (a) - NIR vs. optical stretch parameter. The solid line is the 
best-line fit. The four encircled points are the SNe shown in Fig. Ill (b)- 
histogram of the difference between time of B band maximum estimated 
from the NIR and optical LC fits from the GOLD sample. 


deviated by more than 3cr were removed from the fits. The same 
was repeated for the later data and the two fits were stitched to¬ 
gether. Because our NIR X-correction are reliable up to phases 
—1-40 - 50 days, the data past h- 50 days was not used. We no¬ 
ticed that the NIR stretch ^nir of SNe with iopt ~ 1 was not 1. 
For this reason we additionally stretched the template time axis 
so that the SNe with iopt ~ 1 also have on average ^nir ~ 1- 

Armed with the new templates the above steps were repeated 
to derive our final J and //-band LC templates shown in Fig. 
The templates were additionally normalized to have zero magni¬ 
tude at the first NIR LC peak so that the fitted magnitude offsets 
are the rest-frame NIR peak magnitudes. The weighted rms scat¬ 
ter of the points around template^ cr^^, is cTv/ - 0.02 mag and 
cr^ ^ 0.03 mags for J and //, respectively. The small scatter 
shows that the templates describe the LCs of the GOLD sample 
very well. Between -I-IO and -1-50 days the scatter is considerably 
larger, ~0.13 mag and 0.09 mag for J and H, respectively. 

The GOLD SN sample was fitted with the new templates. 
The NIR LC parameters, together with those derived from the 
optical photometry are shown in Table Figure shows the 
NIR vs. the optical stretch parameters. The best weighted least- 
squares linear fit to the points is shown with the solid line and it 


® weighted by their photometric errors 
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Fig. 7. Example of NIR LC fit of SNe from the GOLD sample. Solid 
line - best fit; dotted lines - fits with = iop,. J band is shown with 
blue and H with red color. 


essentially coincides with the one-to-one correspondence ^nir = 

■^opt- 

Careful inspection of the fits, however, revealed that the scat¬ 
ter in Fig. 1^ arises from real differences between optical and 
NIR stretch. In Fig. [7] are shown the fits of four SNe. The solid 
lines are the fits with ^nir as a free parameter and the dotted 
lines with ^nir = ■Sopt- It is evident that the LCs of SN2002bo 
and SN 2006et are better fitted with snir as a free parameter. 

The zero time of the NIR LC templates is set at Thus, 
the NIR LC fits with the stretch and time of "maximum" as free 
parameters provide an independent estimate of . In Fig. 
is shown the histogram of the difference between the es¬ 

timates of from the NIR and the optical LC fits. The distri¬ 
bution peaks close to zero and has width around 0.5 days. The 
intrinsic scatter of CTt, was estimated by requiring that;^^^ 

with respect to the weighted mean of AfB„„: 




z 





(4) 


was equal to Nsn - 1- Here cr^ is the uncertainty of Afg , 

and Ats^ is computed with weights inversely proportional to 
the denominator in Eq.|^ This way cr, = 0.44 days was obtained 
and the weighted mean of Atg^ was close to zero. 


3 . 5 . Fitting under-sampled NIR LCs 

The analysis of the GOLD sample further reinforced the results 
of Krisciunas et al.] ( 2004c|bl l and |Wood-Vasey et al. ( |2008 1 who 
showed that smooth LC templates stretched according to the op¬ 
tical stretch accurately describe the rest-frame J and //-band 
light curves within 10 days from B-band maximum. Thus, if iopt 
and are estimated from optical observations, by fixing the 
zero time of the NIR template to and setting snir = Sopt, one 
could in principle measure the SN la NIR peak magnitudes with 
a single observation. We used this approach to estimate the peak 
magnitudes of the SNe in our sample and all published SNe with 
under-sampled LCs with observations before H-10 days. This also 
includes the|BN12|sample. However, note that most SNe in this 
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sample have several observations around maximum and the LC 
quality is closer to that of the GOLD sample; often the only dif¬ 
ference is in the photometry uncertainties, which are smaller in 
the GOLD sample. The parameters derived from the LC fits are 
given in Tables [9 -[^ 

When only a single NIR observation is available, depending 
on its phase, possible time offset between the optical and the NIR 
LCs and the difference between snir and Sopt (Fig- will intro¬ 
duce additional uncertainty in the estimated peak magnitudes. 
To estimate this uncertainty a series of simple Monte Carlo sim¬ 
ulations were performed. To do that it was necessary to have an 
estimate of the intrinsic scatter ctj of ^nir around the line defined 
by the snjr = Sopt relation. The scatter cr^ = 0.05 was estimated 
in a similar manner to the estimation of cr, in the previous section 
and the offset from the snir = -Sopt line was essentially zero. For 
given stretch sq 10000 random pairs of stretch and time offset 
(■so.i'j to,i) were generated with the relations sqj = io + N{Q, cr^) 
and toj = N(Q, cr]), where Niji, cr^) is Gaussian distribution with 
mean p and variance cr^. The J and H templates were stretched 
and shifted with each pair (sq,/, fo,/) and the template magnitude 
at phases from -7 to -H10 days were estimated. Then the standard 
deviation with respect to the magnitudes of the unperturbed tem¬ 
plates was computed for each phase. The results for a SN with 
■Snir = 1.0 at several phases are shown in TablefT^ When the ob¬ 
servation is obtained with a few days from the NIR maximum the 
peak magnitudes are fairly insensitive to the uncertainties of the 
stretch and the time of maximum with the scatter being < 0.024 
mag. In the H band the scatter is also small for phases up to -H10 
days. However, in the J band the scatter increases considerably 
beyond phase 0 days reaching 0.12 mag at -I-IO days. Although 
these results were obtained for LCs with only one point, they 
should also be valid for LCs with 2-3 points, provided the obser¬ 
vations are clusters around given phase, e.g. if they are obtained 
with cadence of one day. 

Another set of more realistic simulations based on the LCs 
of the GOLD SN sample were also performed. The peak magni¬ 
tudes estimated from the full LC fits were assumed to be the 
"true" SN peak magnitudes. Then two points were randomly 
drawn from each LC , one from the points within 3 days from the 
NIR maximum, and another one from the points taken later than 
phase h- 4 days. The selected points were fitted with the method 
for under-sampled LCs, and the weighted mean and RMS of the 
difference from the "true" peak magnitudes was computed. The 
mean values were in all cases close to zero indicating that the LC 
under-sampling does not lead to bias of the estimated peak mag¬ 
nitudes. The weighted scatter for the points close to maximum 
was ~ 0.025 mag in both bands. The scatter increases consider¬ 
ably with the points obtained after phase -i-4 days; 0.11 mag in J 
and and 0.07 mag in H. 

Based on the results of the simulations the literature sample 
with under-sampled LCs was split into to groups (i) SILVER 
with the SNe that have observations within 2-3 days from the 
expected time of the NIR maximum light, and (ii) BRONZE with 
the remaining objects. Table [T3] shows statistics of the NIR SNe 
samples. 

Finally, we describe the procedure used to propagate the un¬ 
certainties of the stretch and time of maximum to the uncer¬ 
tainty of the peak magnitudes. We used the following Monte 
Carlo procedure. For each SN, 10000 random pairs of stretch 
and time of maximum {s, t) were generated according to the re¬ 
lations i = iopt +N(0, CTj -i-Bigp,) and t - +N{Q, cr]-i-dtg ), 

where iopt, d-Sopt and are the optical LC stretch and 

time of B-band maximum and their uncertainties, and cr, = 0.44 
days and ctj = 0.05 as determined previously. The LC was then 






















Table 8 . Light curve parameters of the GOLD sub-sample 
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Notes. see text for the references for the NIR photometry. *** The uncertainties are in units of 0.001 mag. The two numbers are the uncertainties of the peak magnitude and the Al-corrections. 
estimated from our ^opt - AM 15 relation. 
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Table 9. Supernova fit parameters. The magnitude uncertainties are given in units of 0.001 mag. 


SN 

SN 2007hz 
SN 2007ih 
SN 2007ik 
SNF20080510-005 
SNF20080512-010 
SNF20080516-000 
SNF20080516-022 
SNF20080517-010 
SNF20080522-000 
SNF20080522-001 
SNF20080522-004 
SNF20080522-011 
SNF20080606-012 
SNF20080610-003 


J c,d 


20.274 (71/49) 
21.375 (101/65) 
21.410(232/47) 
19.742(122/34) 
18.706 (130/28) 
18.865 (67/31) 
19.150 (58/31) 
20.348 (125/42) 
17.684 (33/21) 
17.926 (132/23) 
19.286 (237/40) 
17.284 (25/17) 
18.983 (57/33) 
19.794 (77/37) 


20.679 (125/80) 
21.443 (190/100) 

19.569(115/45) 
19.014 (97/35) 
19.101 (71/39) 
19.214(90/40) 
20.907 (277/68) 
17.887 (30/26) 
18.119(100/28) 
19.756 (337/58) 
17.521 (43/22) 
19.369 (91/42) 
19.811 (92/51) 


gmax^ V at 


350.86 (0.25) 
352.61 (0.42) 

356.78 (0.28) 

604.39 (0.25) 
602.88 (0.35) 
606.33 (0.51) 

609.79 (0.31) 
605.97 (0.81) 
621.17(0.13) 
602.83 (0.45) 
606.33 (2.88) 

616.40 (0.15) 
631.93 (0.50) 
633.43 (0.75) 


_ ■^opt _ 

1.069 (0.020) 
1.184 (0.046) 
1.090 (0.030) 
1.177 (0.029) 
0.743 (0.012) 
1.091 (0.015) 
0.947 (0.015) 
1.003 (0.029) 
1.092 (0.007) 
1.138 (0.021) 
1.260 (0.067) 
1.052 (0.009) 
0.954 (0.017) 
1.047 (0.019) 


AMis" 

0.919 (0.026) 
0.777 (0.046) 
0.903 (0.038) 
0.782 (0.033) 
1.677 (0.047) 
0.902 (0.025) 
1.150 (0.030) 
1.027 (0.069) 
0.901 (0.021) 
0.832 (0.039) 
0.737 (0.040) 
0.935 (0.023) 
1.138 (0.036) 
0.941 (0.032) 


19.555 (18) 
20.767 (25) 
20.617 (27) 
18.826 (22) 
18.038 (49) 
18.165 (38) 
18.379 (21) 
19.692 (52) 
17.027 (07) 
17.422 (49) 
18.991 (185) 
16.653 (09) 
18.684(36) 
18.865 (51) 


19.640(14) 
20.752 (27) 
20.662 (28) 
18.860(18) 
18.114(30) 
18.275 (22) 
18.431 (13) 
19.625 (28) 
17.023 (06) 
17.456 (27) 
19.317(100) 
16.724 (09) 
18.736 (22) 
18.957 (30) 


Notes. JD-2454000. ® estimated from our .Sopt - AMjs relation. The uncertainties are in units of 0.001 mag. The two numbers are the 
uncertainties of the peak magnitude and the k'-corrections. 


Table 11. SN fit parameters for the BN12 sample. 


SN_ ZCMB ./max^’^^ V at 


PTF09dlc 

0.06656 

PTFlOhdv 

0.05424 

PTFlOhmv 

0.03272 

PTFlOmwb 

0.03117 

PTFlOndc 

0.08168 

PTFlOnlg 

0.05592 

PTFlOqyx 

0.06515 

PTFlOtce 

0.03971 

PTFlOufj 

0.07613 

PTFlOwnm 

0.06440 

PTFlOwof 

0.05131 

PTFlOxyt 

0.04831 


18.663 (28/29) 
18.326 (68/24) 

17.137 (79/14) 
19.283 (47/34) 
18.674 (36/25) 
19.028 (40/29) 
17.871 (25/19) 
19.266 (61/32) 
18.731 (30/29) 
18.391 (32/24) 
18.406 (37/23) 


18.948 (26/37) 
18.520 (51/30) 
17.507 (22/20) 
17.401 (46/20) 
19.430 (52/44) 
18.656 (22/31) 
19.200 (35/36) 
18.016 (20/24) 
19.354 (40/41) 
18.997 (29/36) 
18.581 (30/30) 
18.416 (86/29) 


073.79 (0.14) 
342.41 (0.82) 
351.69 (0.03) 
390.32 (0.11) 
389.99 (0.16) 
391.43 (0.27) 
426.36 (0.15) 
441.92 (0.27) 
456.89 (0.13) 
476.68 (0.08) 
473.81 (0.22) 
490.85 (0.11) 


0.991 (0.019) 
1.125 (0.035) 
1.130 (0.002) 
0.885 (0.010) 
1.091 (0.010) 
0.968 (0.020) 
0.913 (0.012) 
1.088 (0.012) 
1.070 (0.021) 
1.090 (0.008) 
1.038 (0.011) 
1.117(0.021) 


1.056 (0.050) 
0.855 (0.056) 
0.846 (0.021) 
1.246 (0.033) 
0.902 (0.022) 
1.110(0.047) 
1.191 (0.026) 
0.905 (0.023) 
0.919 (0.027) 
0.903 (0.021) 
0.954 (0.027) 
0.868 (0.038) 


17.911 (52) 18.051 (21) 
17.492(21) 17.552(22) 
17.278(46) 17.116(24) 
16.709(32) 16.727(11) 
18.427 (8) 18.487 (6) 

18.561 (26) 18.445 (19) 
18.213(16) 18.318(11) 
17.177(14) 17.096(11) 
18.404(48) 18.529(21) 
18.152(9) 18.129(6) 

17.943 (35) 17.826 (08) 
18.452(24) 18.228 (24) 


Notes. JD-2455000. ® estimated from our iopt - AM 15 relation. The uncertainties are in units of 0.001 mag. The two numbers are the 
uncertainties of the peak magnitude and the A'-corrections. 


Table 12. Peak magnitude uncertainties of a SN with %ir = 1.0 when 
only single observation is available. The uncertainties of the stretch and 
the time of maximum were assumed to be 0.05 and 0.44 days, respec- 


Table 13. Statistics of the NIR SN samples used in this work: number 
of SNe, mean redshift, standard deviation around the mean and 

the median redshift. 






Sample 

A^LCs 

< 2cmb > 


median 

Band 


phase of observation 


GOLD 

45 

0.0148 

0.0095 

0.0147 


-4 

0 4 7 

10 

SILVER 

16 

0.0232 

0.0113 

0.0239 

J 

0.013 

0.024 0.049 0.090 

0.123 

BRONZE 

29 

0.0284 

0.0165 

0.0266 

H 

0.006 

0.020 0.008 0.020 

0.006 

GOLD (z > 0.008) 

30 

0.0196 

0.0081 

0.0173 





SILVER (z > 0.008) 

15 

0.0244 

0.0106 

0.0240 





BRONZE (z > 0.008) 

28 

0.0293 

0.0162 

0.0268 

fitted with each perturbed pair (i, f) and the standard deviation of 

BN12 

12 

0.0556 

0.0161 

0.0559 

the 10000 peak magnitude estimates was added 

in quadrature to 

Our new sample 

14 

0.0942 

0.0448 

0.0850 

the uncertainty of the peak magnitude derived by the LC fitter. 






4. Results 




function of AM 15 with 

a slope of ~ 0.11. 

The parameters of 


4.1. Colors at maximum light 

4.1.1. (B-y)B_ 


the (B - V)b^,q - AMi 5 relation were determined by linear fit 
of the bluest objects (see also Folatelli et al.pOl^l. The disad- 


Phillips et al. (19991 were the fist to suggest that the intrinsic 
B - V color of SNe la at maximum {B - V)b„„„ is a linear 


vantage of this approach is that it requires a prior knowledge of 
which SNe are likely to be unreddened and that only a small 
fraction of all available data is used. Here we followed a differ¬ 
ent path to study the distribution of observed {B - V)b„.^ colors 


Article number, page 12 of 38 


























Stanishev et al.: Type la Supernova Cosmology in the Near-Infrared 


with the much larger sample available today. [Nobili & Goobar| 
(|2008]l have shown that the intrinsic SN la {B - V)b^„ color 
has a dispersion of about 0.05 mag. The observed colors may 
be further reddened by interstellar dust along the line of sight, 
evidence for which is clearly seen as an extended red tail of the 
color distribution in the nearby SN la sample. Let us assume as 
in|Jha et al.|(|2007|l that (i) the intrinsic {B - V)b^^ colors follow 
normal distribution with mean cq and standard deviation ctc, and 
(ii) that the reddening E{B - V) follows exponential distribution 
with scale for E{B - V) > 0 mag and is zero otherwise. Un¬ 
der these assumptions the probability density function (PDF) of 
the observed {B - V)b^„ color, p({B - y)^„.„;co Q-cTe), is the 
convolution of the two functions (|Jha et al.|2007t |Marriner et al.| 

|20TlT l: 


p(c;Co,Crc,Te) oc 


f 


exp - 


(c - Co - t)^ 


2(Tr 


exp 


-—\dT, 


(5) 


where c is the observed {B - V)b^ color and the integration is 
over the color excess from extinction, t. 

Figure IM shows the (B - colors measured with uncer¬ 

tainty less man 0.04 mag from 393 LCs. The dotted line has a 
slope of 0.12. It describes the blue edge of the distribution of the 
SNe with AM 15 <1.4 well, but beyond that the colors are clearly 
increasingly redder than the line. It is clear that a non-linear 
function is need to describe the average intrinsic {B-V)b^^ color 
as a function of AM 15 . We tried different functions and found that 
a simple exponential function worked best, while a second-order 
polynomial was not flexible enough to describe the data. 

The parameters of the exponential function that describes 
{B - were determined by fitting Eq.|^to the histogram of 

E{B - y)obs = (B - - (B - V)b^,o computed in 0.05 mag 

wide bins. Unfortunately, there is no unique way to define the 
quantity to be minimized. We tried several approaches, which 
are appropriate under the reasonable assumption that cTc is inde¬ 
pendent on AM 15 . First the histogram of the whole data set was 
fitted and the parameters that give the lowest cr^ was selected. 
The reasoning behind this is that with incorrect parameters the 
histogram will be additionally blurred and hence a larger cr^ will 
be derived; the correct parameters should thus give the lowest 
(Tc. However, this approach was unsuccessful and the blue edge 
of the distribution of the SNe with AM 15 >1.4 was not described 
well. The likely reason for this is that the number of SNe with 
AMi 5 > 1.4 is small and the final solution is mostly driven by 
the SNe with AM 15 < 1.4. The approach that we found to work 
best to address this problem was to split the sample into three 
bins - SNe with AM 15 < 1.1 (182 points), 1.1 < AM 15 <1.5 
(163 points) and AM 15 >1.5 (48 points) - and fit histogram of 
E(B — y)obs in each bin independently. In the fitting the co pa¬ 
rameter in Eq.j^was treated as an unknown constant. It is easy to 
see that with the correct {B - V)b„„,o - AM 15 relation ctc should 
be the same in all three bins and the cq’s should be close to zero. 
Thus, the quantity to minimize was the scatter of ctc around the 
mean in the three bins and the scatter of cq around zero added in 
quadrature. The Te parameters were not required to be the same 
in the three bins because the distribution of the reddening may 
be different. Eor example, the SNe with larger AM 15 are fainter 
and this may lead to bias against discovering reddened SNe. 

With the above procedure the following function was de¬ 
rived; 


(B - V)b„„,o = 0.036 exp 


AMi 5 - 1.3 
0.3455 


0.084(+0.01). ( 6 ) 


This relation is shown with the solid line in Eig. IM and the his¬ 
tograms in the three bins and the best fits in EigT]^. The mean 



(B-V)b -(B-V)b 

' ' “max ' ' “max,( 


Fig. 8 . (a) The ( 6 - color of the SNe used in this work as function 
of AM 15 . The solid line shows the relation defined by Eq.l^ The dashed 
line has a slope of 0.12 and is shown for a comparison; (B) histograms 
of (S - y)B„„ - (fi - y)B„„,o for SNe with AM 15 < 1.1 (182 objects), 
1.1 < AMi 5 < 1.3 (163 objects) and AM 15 > 1.5 (48 objects) with the 
correspondign fits with Eq.|^ 


value of cTc in the three bins is 0.045 mag with estimated un¬ 
certainty of 0.01 mag. The uncertainty of the free term in Eq.|^ 
+ 0.01 mag, was estimated from the fits of the histograms, which 
give uncertainty of ~ 0.01 mag of cq. The constant 1.3 was de¬ 
rived empirically. Several values between 1.1 and 1.4 were tried 
and the one providing the solution with smallest cr^ was selected. 
Regarding Tg, the following values were derived in the three bins, 
Te = 0.14 + 0.02 mag, Te = 0.13 + 0.02 mag and Tg = 0.11 + 0.04 
mag. We note that the results are somehow dependent on the bin¬ 
ning on AMi 5 and the color, and color error cut (0.04 mag in our 
case). When a larger sample of SNe, especially with AM 15 > 1.4, 
becomes available it will be possible to estimate the parameters 
more accurately. 


The values of o-g and Tg that we find are close to those of 
|Jha et al.| ( |2007[ ), [(janeshalingam et al.| ( |2010 l and Eolatelli et al.| 
( |2010| l, which were estimated from the late-time SN la colors at 
+35 days using the Lira-relation ( |Lira| 1995] [Phillips et al.|199^ . 
The Lira-relation is considered as one of the best methods to es¬ 
timate the host galaxy reddening of SNe la. In another analysis 
of late-time B - V color curves of SNe la [Burns et al.j ( [2014| l 
introduced a new stretch parameter of the B - V curves, sbv- 
The authors estimated the intrinsic scatter of the B - V color to 
be 0.06 - 0.09 mag. The intrinsic scatter of (B - V)b^ that we 
estimate is smaller than the scatter of B - V at late times de¬ 
rived by the above authors. This indicates that at least for normal 
SNe la the (B - V)b„,^ color should provide as good estimate 
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Fig. 9. (a) Jmax-fimax color VS. AM 15 . The solid shows the least-squares 
linear fit to the GOLD sample; (b) vs. E(B - V')obs- The solid 

line is the least-squares linear fit to the GOLD sample, the dashed and 
the dotted lines are the expected relation for Milky Way-type dust and 
Rb = 2.7 and 4.1, respectively. 


of the reddening as the Lira-relation. This may be important for 
high-redshift studies where the observations do not extend to late 
phases and the extinction must be estimated from the colors at 
maximum. 

As a final note, we would like to point out that the width of 
the observed color distribution is not determined only by the in¬ 
trinsic color scatter cTc and and the distribution of the reddening 
(here parametrized by Te). Several other other factors, such as 
the observational errors, the fact that S -corrections were not ap¬ 
plied to all data and possible deviation of the fitting LC templates 
from the real SN LCs will all increase the width of the observed 
distribution and lead to overestimate of cTc and Tg. Some of these 
sources of errors may be different for the different samples and 
therefore ctc should be regarded as a mean value for the whole 
nearby SN la sample. The mean uncertainty of (B - in the 
sample used in this analysis is ~ 0.015 mag and taking it into 
account the true intrinsic scatter ctc may be ~ 0.04 mag. 


4.2. Jmax - Hmax, J max and Vb„^ - H^ax 

Figure|^shows the J^^x — ^fmax color vs. AM 15 and E{B - VOobs- 
The 7max - H^ax color shows a fair amount scatter, largely due to 
our and the BRONZE samples. The scatter in the GOLD sample 
is rather small and the weighted least-squares linear fit indicates 
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a possible weak dependence on AM 15 with a slope of 0.14 +0.06. 
The fit of the GOLD sample implies that a typical SN la with 
AMi 5 = 1.1 has Jmax “ Hmax - “0.23 + 0.02 mag. ¥or !DoF 
of the fit to be one it is necessary to add 0.08 mag in quadrature to 
the errors, which can be interpreted as the intrinsic scatter of the 
Jmax - Hmax color. Similar intrinsic scatter is also inferred from 
the BN 12 sample, however, /max - Hmax color in this sample is 
systematically redder by about 0.1 mag. 

The GOLD sample also reveals a weak linear relation be¬ 
tween Jmax - Hmax and E(B - y)obs- The relation is shown with 
a solid line in Fig. [^. The expected relations for Milky Way- 
type dust with Rv - 1.7 and 3.1 are shown with the dashed 
and the dotted lines, respectively. The slope of the best-fit line, 
0.12 + 0.03, is consistent with a low Ry value around 1.5. The 
fit indicates that an un-reddened SN should have Jmax - Hmax = 
-0.25 ± 0.03. We also note that there is no apparent dependence 
of Jmax - Hmax of the redshift. For both the GOLD sample with 
mean redshift z - 0.02 and for our sample with mean redshift 
z = 0.09 Jmax - Hmax colof is -0.24 mag. 

Our estimate of the intrinsic Jmax - Hmax color of a SN with 
AMi 5 = 1.1 mag is similar to the results of |Krisciunas et alT 


(|2004c I and |Wood-Vasey et al.| ( |2008) . It is also consistent with 
Burns et aL| ( |2014| l who, using a different approach, found that 
a SN with stretch equal to one has /max - Hmax - -023 mag 
with intrinsic dispersion of 0.09 mag. It is interesting to note, 
however, that |Folatelli et al. (2010 1 based on part of the data we 
used in our analysis find it to be close to zero. From Fig.|^one 
can see that this is inconsistent with our results. We also note 
that the observed variation of the Jmax - Hmax color is far too big 
to be caused by dust extinction of any reasonable law. 

Finally, Fig.[^shows VB^-Jmax and Vb^,,-H max colors vs. 
E(B - y)obs- Most of the SNe with E(B - y)obs ^ 0.4 mag align 
along the line tracing Milky Way-type extinction with Ry = 1.7. 
There is also a set of SNe with E{B - y)obs between 0.1 mag 
and 0.6 mag that seem to follow the line with Ry = 3.1. The 
two SNe with E{B - y)obs < -0.15 mag appear to be outliers. 
However, one of them SNF20080522-004 has large uncertainty 
of the maximum B - V colors of ~ 0.2 mag. The other object 
is SN 2006hx, which is ~ 1 mag outlier on the optical Hubble 
diagram, hence its anomalous y-NIR colors. 


4.3. SNe la as standardizable candles 


In this section we use our estimates of the SN peak magnitudes 
and AMi 5 parameters to re-evaluate the precision of SNe la as 
standardizable candles. The main goal of this paper is to do that 
in the NIR J and H bands. However, because several new large 
sets of optical observations of nearby SNe la have become re¬ 
cently available we also analyzed the B and V bands. 

To standardize the SN la peak magnitudes in band X, mx, we 
adopt the following two-parameter relation; 

m*^^mx-axX (AM 15 - 1.1) - x E(Y - Z), (7) 

where E(Y - Z) is the color excess in bands Y and Z caused by 
dust and which can be converted into absorption in band X via 
R^^ X E(Y - Z). Throughout this work we use E{B - V), which 
was calculated with the {B - y)B„,„,o - AM 15 relation derived in 
the previous section, and omit the superscript index of Rx. 

The ax and Rx parameters, and the absolute magnitude in 
the band X, Mx, of a fiducial SN la with AM 15 = 1.1 were deter¬ 
mined by ;^^^-minimization of; 


A 


A^set A^SN 


'Iv ■ ■ 

Yij 


■M' 


. , . , (T? . . + cr^ . + (T^ . . 
i=l j=l lc,ij intr,; ext ,/,7 


( 8 ) 
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Fig. 10. Vb^,^ -Tmax (left) and Vb^,^ -ffmax (right) colors vs. E(B- F)obs- 
The dashed, the solid and the dotted lines show the expected relations 
for Milky Way-type dust and Rb = 2.7, 2.9 and 4.1, respectively. 


with iterative 3cr outlier rejection. In Eq.j^we allow for the dif¬ 
ferent data sets to have different absolute magnitudes M'^. Here 
(omitting the summation indices) (Text = 150 km s“' accounts for 
the random host galaxy peculiar velocities and 


2 2,22 , d2 2 

O'lc - ‘^mx + “xCTaM,, + 


(9) 


The nearby SNe la sample by its own does not constrain the cos¬ 
mological parameters and the distance modulus fj. was computed 
with the adopted cosmological parameters. 

The terms CTintr,/ are included in order to account for possible 
intrinsic SNe la luminosity scatter, which cannot be corrected 
for by the model defined by Eq. [7] Eollowing |Marriner et al.| 
(|201 l|l we allow for the possibility that (Xintr,/ includes a non¬ 
zero intrinsic color scatter (Xcoi.intr,,, which does not correlate with 
luminosity. Thus 


2 _ 2 , n2 2 

^intr,/ “ ^X^cotintr,/’ 


( 10 ) 


where cr,„x,inti,; is the intrinsic SN la luminosity scatter. In princi¬ 
ple, (Tmx.intr,/ includes also contribution from other sources of un¬ 
certainty such as calibration uncertainties, possible mismatch be¬ 
tween the templates and the actual fitted LCs, etc., which could 
be different in different data sets. Eor this reason, the values of 
the terms (T,nx,mtr,i were obtained so that the reduced is one 
for each data set. This approach helps to prevent samples with 
poorer-quality data from decreasing the accuracy of the derived 
parameters (Amanullah et al.|2010|l. Eor low-redshift samples 


the inferred value of cr, 


wx.intr,/ 


is also dependent on the adopted 


value of the dispersion of the galaxies peculiar velocity (Text, 
which is somehow uncertain. In general, the intrinsic color scat¬ 
ter (Tcoi.intr,/ could be different for each data set, but in our analysis 
we assumed that it was the same for all data sets. 


4.3.1. B-band 


The B-band peak magnitudes of the whole nearby SN la sam¬ 
ple with redshift above 0.015 plus the 16 SNe from this work, 
348 peak magnitudes in total, were fitted with and without in¬ 
trinsic color scatter. The value of intrinsic color scatter was fixed 
to (Tcoi.intr = 0-04 mag, the value which was derived in Sec. 14.1 
from the analysis of the color distribution. The derived 

parameters are listed in Table [l^ and the results of the fit with 


o'coi.intr = 0-04 mag are shown in Eig. 11 The overall weighted 


1(T scatter around the best fit with color scatter is 0.185 + 0.007 
mag. In the individual samples the scatter varies between 0.14 
and 0.30 mag and the corresponding intrinsic luminosity scatter 
cTB.intr N also different. Erom Table [l4| one can also see that the 
absolute magnitudes Mg for the different data sets span a range 
of ~ 0.16 mag, possibly indicating significant systematic differ¬ 
ences between the data sets. The parameters of our sample fall 
well in the ranges of the other SN samples. Note that in this pa¬ 
per as a measure of the scatter of the Hubble diagram we use 
the weighted standard deviation. The weights are inversely pro¬ 
portional to the term in the denominator of Eq. Although we 
quote the observed scatter with no explicit attempt to correct for 
systematic effects, galaxies peculiar velocity term is included in 
the weights and thus the more nearby SNe are given less weight 
in both the Hubble diagram fits and when computing the scatter. 

It is interesting to compare our results with those obtained by 
Betoule et al. ( 2014| l with the popular SALT2 EC fitter and us¬ 
ing similarly heterogeneous low-redshift samplj^ Eor this we 
selected only the 104 SNe that are in common. A direct com¬ 
parison shows that the SALT2 xl and our stretch parameters are 
linearly related as iopt = 1.0 -H 0.108 jd. The comparison of B^ax 
shows only a small offset with a median value of 0.022 mag and 
a robust estimation of the scatter of 0.042 mag. The median dif¬ 
ference between our (B - color and SALT2 c parameter 

is 0.047 mag with a robust scatter of 0.042 mag. Both differ¬ 
ences show small trend with the (B - V)b„,„ color with slope 
0.17 + 0.04 for the B^ax difference and 0.24 + 0.04 for the color 
difference. If we restrict our analysis only to the SNe in common 
with Betoule et al. (2014 1 we find the weighted scatter of the 
Hubble residuals to be 0.15 mag, which compares well with the 
0.14 mag that we estimate from the data published by |BetouTe| 
et al. ( |2014[ ). The intrinsic luminosity scatter cannot be directly 


compared because we assumed it to be different for each sample. 
[Betoule et al.| ( |20f4 l estimate it to be 0.12 mag and we find the 
same value for the CSP and Calan-Tololo samples, and larger in 
the others. Given the significant differences between the two EC 
fitting approaches, we conclude that the results compare rather 
well. Eurther, analyzing our whole sample, but with a color cut 
of (B - y)B„„ < 0.3 mag, we find overall scatter of 0.18 mag. 
Compared to the sample of Betoule et al.| ( |2014[ ), our color-cut 
sample has about 150 new SNe from Berkeley and CfA4 sam¬ 
ples, and about 60 more SNe from CfA3 and CSP samples. Prom 
Table [T4| one can see that Berkeley, CfA3 and CfA4 samples all 
have observed scatter > 0.2 mag, which appears to be the likely 
reason for the larger observed scatter in our analysis. 


The sample used by 


Betoule et al. 


These authors also applied a color cut (B - V)b, 
not include color scatter in their model. 


1 2014 1 is contained in our sample. 

0.3 mag and did 
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Table 14. Results from the analysis of the optical observations. For each sample are listed the absolute B magnitude, the systematic luminosity 
scatter, the weighted observed scatter, and the number of used and total SNe. 


Data set 

Mb 

O" B.intr 

_ b 
^^obs 

-^used(-^tot) 

Mb 

O" B.intr 

*^obs 

^used(^tot) 


Color scatter model CTcoi.intr = 0.04 mag 

Pure luminosity scatter model CTcoi.intr 

= 0.0 mag 

Berkeley 

-19.239 (0.023) 

0.141 

0.201 

87 ( 88 ) 

-19.211 (0.022) 

0.179 

0.196 

87 ( 88 ) 

CSP 

-19.378 (0.026) 

0.095 

0.163 

49 (49) 

-19.326 (0.024) 

0.129 

0.142 

48 (49) 

CfAl 

-19.369 (0.078) 

0.151 

0.232 

9(9) 

-19.336 (0.073) 

0.178 

0.217 

9(9) 

CfA2 

-19.250 (0.052) 

0.130 

0.191 

14(14) 

-19.226 (0.053) 

0.180 

0.195 

14(14) 

CfA3 

-19.325 (0.023) 

0.095 

0.184 

84 (89) 

-19.268 (0.025) 

0.183 

0.212 

87 (89) 

CfA4 

-19.272 (0.041) 

0.185 

0.251 

39 (39) 

-19.248 (0.041) 

0.221 

0.251 

39 (39) 

Kowalski 

-19.147 (0.083) 

0.125 

0.204 

6 ( 6 ) 

-19.149 (0.070) 

0.140 

0.173 

6 ( 6 ) 

Krisciunas 

-19.340 (0.048) 

0.0“ 

0.102 

8 ( 8 ) 

-19.311 (0.031) 

0.060 

0.086 

8 ( 8 ) 

Calan-Tololo 

-19.311 (0.043) 

0.088 

0.181 

18 (18) 

-19.292 (0.041) 

0.137 

0.169 

18 (18) 

BN12 

-19.155 (0.047) 

0.050 

0.162 

12 ( 12 ) 

-19.133 (0.050) 

0.142 

0.171 

12 ( 12 ) 

This work 

-19.214(0.053) 

0.108 

0.210 

16(16) 

-19.211 (0.052) 

0.160 

0.208 

16(16) 

Mean Mb 

-19.291 (0.070) 


0.185" 


-19.261 (0.059) 


0.180" 


Gx 

2.916 (0.064) 




2.692 (0.062) 




Rx 

0.720 (0.038) 




0.698 (0.037) 





Notes. with trg, intr = 0.04 mag the reduced = 0.6. weighted standard deviation. weighted standard deviation of all points around the 
Hubble line. Its uncertainty is 0.007 mag. 



0.01 


0.10 



0.8 1.0 1.2 1.4 1.6 1.8 



Fig. 11. 6 -band Hubble diagram. The right panels show the Hubble residual when only the LC shape or the extinction correction is applied. 


IScolnic et al.| ( |2014| l have claimed that if the color scatter 
is erroneously treated as luminosity scatter this will lead to bi¬ 
ased estimates of ub and Rb coefficients. Rb in particular can be 
biased by as much as 30% toward lower values. These authors 
claim that when the color scatter is treated properly, they recover 
the usual Rb - 4.1 instead of the much lower value found in 
other SN la analyses. In our analysis with and without intrinsic 
color scatter we obtain very similar a^’s. In the analysis with 
zero color scatter we find Rb to be by ~ 8 % lower than in the 
analysis with color scatter of (Tcoi.inti = 0.04 mag (Table [T4|. In 
both cases the value that we recover is Rb ~ 2.7 - 2.9. Repeating 
the analysis only with the SNe with E{B - y)obs < +0.4 does not 
change these values much. From the lower right panel of Fig. [TT| 
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one can see that Rb — 4.1, shown with the dotted line, is incon¬ 
sistent with the data. 


To assess the reliability of the estimation of ob and Rb we 
performed a series of Monte-Carlo simulations. Artificial data 
sets that have the same number of SNe as the full data set, the 
same distributions of {B - V)b„„, AM 15 and redshift were gen¬ 
erated for several values of gb, Rb, luminosity and color scatter. 
The results confirmed the presence of the bias in Rb described 
by Scolnic et al. (20141 when the color scatter is treated as pure 
luminosity scatter. However, the bias was only at the level of 
< 10% and even less in gb- The differences between our results 
and those of Scolnic et al.| ( 2014| l are likely due to the different 
color distributions of the SN samples analyzed. The majority of 
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the SNe in the sample simulated by |Scolnic et al.| ( |2014 1 have 
(B - V)b^^ < +0.2 mag and none with (B - > +0.3 mag, 

while the full nearby sample has many SNe with (B - V)b„.^,_ 
mag in the range 0.2-1.5 mag. The larger the color range the 
better is the constrain on Rb and there is less room for bias. For 
example, in our case 30% bias means that Rb should be 3.9 in¬ 
stead of 2.9, but from the lower right panel of Fig. one can 
see that Rb - A- describes the data rather poorly. It should be 
emphasized, however, that the simulated samples have single Rb 
for all SNe, which most likely is not the case for any real SN 
sample. Indeed, in a recent study of optical and NIR photom¬ 
etry of a nearby SNe la sample Burns et al. ( |2014 l find large 
variations of Ry from event to event. The same conclusion was 
reached by Amanullah et al. ( 2015| l, in a study that also included 
UV measurements, where the differences in extinction laws are 
the largest. 

Folatelli et al. (|2010l analyzed a sub-sample of the SNe la 


used in our work. They find a stronger dependence of SN la peak 
magnitudes on AM 15 with ub - I, while we find a smaller value 
of qb - 0.7. The reason for this difference is unclear, but we 
note that [Folatelli et al.|(|2010| l determined AM 15 with the linear 
relation shown in Fig. which can affect the estimation of 

The analysis of the color indices at maximum light suggested 
that several SNe may be reddened by dust with Rb - 4.1 instead 
with the mean value of 2.9 that we find for the whole set (see 
Fig. [T0| l. The mean color excess of these SNe is ~ +0.3 mag. 
If these SNe had to be corrected with Rb - 4.1 but we used 
Rb - 2.9 instead, then their Hubble residuals should have been 
on average ~ +0.3 mag. However, this is not the case and the 
Hubble residuals are on average close to zero. Thus, most likely 
the behavior of the maximum colors of these SNe is an intrinsic 
property and not caused by dust extinction with Rb - 4.1. We 
also note that all SNe along the Rb - 4.1 line, except for one. 


belong to the normal-velocity group of SNe la (Blondin et al. 
2012| Wang et al.|2009a 1 . For the whole NIR sample, about 20% 


belong to the high-velocity group and thus the SNe along the 
Rb - 2.9 line are mixture of both velocity groups. This suggests 
that the behavior of the maximum colors is most likely unrelated 
to the velocity dichotomy noted by [Wang et al.|([2009a|l. 


4.3.2. NIR Hubble diagram 

The model defined by Eqs. |7p0| was also used to analyze the 
NIR J and H peak magnitudes of the SNe with z > 0.008. The 
data was considered consisting of five SN samples; the new SN 
observations presented in this paper, the |BN12| sample, and the 
published SN sample, which is split into three sub-samples - 
GOLD, SILVER and BRONZE. Two of the SNe from our new 
sample, SN 2008bz and SNE20080512-008, were observed in 
the NIR at stretch-corrected phase —1-13 days. This is beyond 
the range where our fitting method could be used and these two 
SNe were excluded from the analysis. Further, two SNe from the 
CSP sample, SNe 2005ku and 2007hx, were found to be signifi¬ 
cantly brighter than the best-fit Hubble line in H band. We con¬ 
tacted CSP and they confirmed that most of H band photometry 
of these two SNe were unreliable (C.Contreras, private commu¬ 
nication). Only one H observation of SN 2005ku was deemed 
good and the updated photometry was even brighter than the 
published one. For this reason SNe 2005ku and 2007hx were 
also excluded. The results of the analysis are shown in Table [TS] 
In addition to the fitted values and uncertainties of ax, Rx and the 
absolute magnitudes of the different samples, we also give the 
weighted tms of the whole sample and the five sub-samples, as 


well as the systematic luminosity scatter that makes I do f - 1 

for each of the sub-samples. 

The data were first analyzed assuming that SNe la are natural 
standard candles in the NIR and no corrections to the peak mag¬ 
nitudes are needed. The weighted scatter around the Hubble line 
was 0.155 mag in J and 0.140 mag in H, which is better than the 
scatter we found in our analysis of the optical data. Inspection of 
the Hubble residuals revealed a weak dependence on AM 15 and 
the color excess E(B- y)obs- When the ajjn and By/n correction 
terms were included in the analysis the scatter of the Hubble 
residuals decreased to 0.147 mag in J and 0.135 mag in H. In 
general, the ajjH coefficients are fairly small, which combined 
with the lack of many of highly-reddened SNe, is the reason for 
the rather small improvement. 

Further analysis of of Hubble residuals hinted at possible de¬ 
pendence on the 7max“7/niax color {V—lJjH} were also tested, but 
the highest degree of correlation was found with the /max - T/max 
color). The term: 

-rx X [(y^ax - 7/max) - (J - H\] , ( 11 ) 

was added to Eq. where (7 - H)q was set to the mean 
■/max - //max color of -0.24 mag. Eigure [T^ shows the //-band 
Hubble diagram with the an, Rh and Jh correction terms ap¬ 
plied. Only the //-band Hubble diagram is shown because we 
found that the J- and //-band residuals are highly correlated after 
applying all three corrections. Eigures [T^c-e) show the Hubble 
residuals when only two of the an, Rh and jh terms are applied 
vs. the quantity multiplied by the third, e.g. Eig.fT^e) shows the 
Hubble residuals vs. /max - //max when only the an and Rn terms 
are applied. Figure [T^f) is the same as Fig. [T2|e) but for the J 
band. With all three corrections included the weighted scatter be¬ 
comes 0.14 + 0.01 mag in both bands. This is somehow surpris¬ 
ing because several previous studies have claimed that //-band 
Hubble diagram has considerably smaller scatter than J (e.g., 
jWood-Vasey et aL]|2008| ), but we find that SNe la are equally 
good standard candles in both bands. 

The presence of non-zero y term is interesting because it has 
not been previously reported. However, the inclusion of the y 
term does not lead to a significant decrease of the Hubble resid¬ 
uals. The quantity that is most interesting for cosmology is the 
intrinsic scatter (Xintr, which should be independent on the pho¬ 
tometric uncertainties, as opposed to the scatter of the Hubble 
residuals. From Table [TS] one can see that including the y term 
does not lead to decrease of the inferred crimr in all samples. In J 
band BN 12 and our samples are strongly affected, SILVER and 
BRONZE only marginally, and the GOLD is unaffected at all. 
In H band only SILVER and BRONZE are affected, but the rest 
are not; BN12 and our samples even got worse. It is evident that 
the non-zero y term is not driven by the best observed GOLD 
sample, which casts doubt of whether it is real or an artifact of 
the small sample sizes. However, the range of /max - //max color 
in the best-observed GOLD and BN12 samples is smaller than 
in the others. Thus, the question of whether the Hubble residu¬ 
als depend on the /max - //max color should be answered in the 
future when a larger sample of well-observed SNe la in the NIR 
becomes available. 

The value of (Timr is a very important quantity because it tells 
us to what accuracy SNe la peak magnitudes could be standard¬ 
ized. It is however a quantity that is difficult to estimate as all 
sources of uncertainties of the peak magnitudes have to be prop¬ 
erly accounted for. If that is the case, then all samples should 
give similar estimates of crimi. However, the exact value of the 
random component of galaxy velocities cText is rather uncertain. 
This uncertainty affects the more nearby GOLD, SILVER and 


Article number, page 17 of 


38 





























A&A proofs: manuscript no. nir_sne_submitted 


W) 


O 

D 

'c 


O 


b 


s 

o 

J= X) 

W) 03 


-D 

’bb 


(U 

x: 'Hi 

o 

«t 

03 C 

•V 

^ t: 


’O 

X5 


03 


O' .S 


^ •- 
S-H c/5 

■S C 
‘S 
S 

"S 

^ § 



h 


c 


b 

QJ 


a. 


S 

t/3 


N 

{/I 


U 


s 


o 



c 


b 

(N 


z 

a 

Ui 

m 





c 


b 

w 


N 


z 

B 

o 

Ui 

a 


m 



S 


h 


c 


b 

Pi 


w 

t/3 

> 

a 



on 







c 


b 

Q 



B 

O 

Ui 

O 




-n 


OJ 



t/3 

a 


u 




X 


«s 








Article number, page 18 of 


38 





























Stanishev et al.: Type la Supernova Cosmology in the Near-Infrared 



0.01 0.10 

^CMB 



Fig. 12. //-band Hubble diagram for the full sample with the a^, and ju correction terms included. Panels (c-e) show the Hubble residuals 
when only two of the ajj, and ju terms are applied vs. the quantity multiplied by the third, e.g. in panel (e) the Hubble residuals vs. Jmax “ ffmax 
when only the and Rfj terms are applied. Panel (f) is the same as panel (e) but for the J band. In panel (d) the dashed lines are the extinction 
corrections for Milky Way extinction law with Ry = 1.7. The two red point surrounded by black circles are SN 2008bz and SNF20080512-008, 
which have been excluded from the fit (see text). 


BRONZE samples and will have relatively small effect on BN12 
and our samples. From Table [T5]one can see that the BN12 sam¬ 
ple suggest that cTintr is likely less than 0.1 mag. In J band the 
BN12 sample has very small crma - 0.028 mag (this further dis¬ 
cussed later). In J band our new sample also has rather small 
c^inti- - 0.08 mag but in H is it considerably larger ~ 0.14 mag. 
Interestingly, even though we assumed rather conservative esti¬ 
mate of (Text = 150 km s“', the GOLD sample also has (Tintr - 0.1 
and small observed scatter. 

Although (Tinti from the GOLD, SILVER and BRONZE sam¬ 
ples depends on the assumed value of (Text, the inferred values 
of (Tint! can be compared. In general, (Tintr in the SILVER and 
BRONZE samples is larger than in the GOLD sample by a factor 
of 1.2-1.5. Our sample also has larger (Tintr- Thus, it appear likely 
that in spite of our efforts to identify and include all uncertainties 
that affect the peak magnitudes estimated from incomplete LCs, 
some additional scatter still creeps in. One possible source of 
such scatter is systematic deviation from the knir = iopt relation. 
Careful inspection of Fig. §1 reveals that the points in the range 
■Sopt = 0.92 - 1.02 are systematically above the line, while those 
the ranges Sopt = 0.8 - 0.9 and Sopt = 1.02 - 1.15 are system¬ 
atically below it. This raises the interesting possibility that knir 
has a bi-modal distribution with mean values of knir - 0.8 and 
■snir - 1-05 for the SNe with iopt below and above ~ 0.91, re¬ 


spectively, but clearly more observations are needed to confirm 
whether these systematic trends are real or not. 

Despite the differences of the inferred values of cTi„tr and the 
observed scatter, the SN la absolute magnitudes estimated from 
the GOLD, SILVER, BRONZE and our new samples are in a 
good agreement. This indicates that using under-sampled LCs 
does not lead to significant bias in estimating the NIR peak mag¬ 
nitudes. However, the |Barone-Nugent et al.| ( |2012| l sample has a 
systematic offset and is on average 0.15-0.20 mag fainter. This 
difference is not negligible and appears to be too large to be ex¬ 
plained with differences between the photometric systems used. 


The dependence of the peak magnitudes on AM 15 and E(B - 
V)obs is rather weak. The ax coefficient is ~ 0.16 - 0.20 in both 
bands and is smaller than the estimation of |Kattner et al.|p012| l 
assuming bi-modal distribution of SN la NIR peak magnitudes. 
The method used by Kattner et al. ( 2012| l to estimate AM 15 is the 
same as in |Folatelli et al. 12010 1 , which likely contributes to the 
differences in ax coefficients. In both J and H, the ax coefficient 
is more than twice smaller than in B, which further supports the 
conclusion that SNe la are better natural standard candles in the 
NIR than in the optical. The values of Rjjn that we find are in the 
range 0.2-0.4, which is is very close to the expected values for 
Milky Way type of extinction with low Ry - 1.5-2. In Fig.[T^ 
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the dashed line shows the extinction corrections for Milky Way 
extinction law with Ry - 1.7. 

4.4. V - J and V - H color indexes evolution 

[Krisciunas et al.| ( |2004b| l have studied the V-NIR colors of 
SNe la and found that their evolution between phases -10 and 
+ 10 days is fairly uniform. They also showed that SNe with 
AMi 5 < 1.1 mag and AM 15 > 1.1 follow different V — H loci. 
We studied the evolution of the V - J and V - H colors with 
the SNe of the GOLD sample with E{B - V) < 0.2 mag. The 
observed V, J and H magnitudes, corrected for the Milky Way 
dust extinction, were converted to rest-fram^^and the observed 
epochs were converted to stretch-corrected phases using Eq. 
with Sopt. The V magnitudes were interpolated at the phases of 
the NIR observations. Most SNe have well-sampled V band LCs 
and the interpolation was done with smoothing splines. Only in 
a few cases interpolation from the templates fits was used. We 
assumed that the evolution of the colors indexes with the stretch- 
corrected SN phase is the same for all SNe and that the offsets 
between color curves of the different SNe arise from different 
dust extinction in the host galaxies and possibly from the intrin¬ 
sic color scatter. Because the slope of the "luminosity-light curve 
width" relation in V and the NIR bands is different, it is to be ex¬ 
pected that the offset may also be a function of the decline rate 
parameter AM 15 . Under these assumptions and in the ideal case 
of no intrinsic color scatter, the observed color curves V -X can 
be described by: 

y - X = f^v-x)W + R(V-X) E(B -V) + fl(y_x)(AMi 5 - 1.1), (12) 



Fig. 13. V - J and V - H un-reddened loci for SNe with AMis = 1.1. 
The 3rd-order polynomial fit of the data between phases -11 and 12 
days is also shown. For clarity the error bars are not plotted. 


where X is either J or H, is the average color evolution 

as function of the phase, R(v-X) - Ry — Rx and a[v-x) — ay- ax- 
The function fy-x)(<f>) was assumed to be 3rd-order polynomial. 
To derive its coefficients, and R^y-x) and a(v-x) parameters, the 
following statistic was minimized: 

r ZjZj ^2 +2 • 

i=i j=i ^ (y-x),y inti- 

where 

(V-Xy = iV-X)-Ry^x)E{B-V)-ay^x)(AMi5-l-l) (14) 

and the summation is oven all SNe and observations obtained 
between phases -11 and +12 days. The term cTintr was added so 
that the reduced of the fit was one. 

After the first iteration we noted that 4 SNe had color curves 
significantly deviating from the rest. SNe 2002dj, 20061e, 20061f 
and 2008gp were all too blue at all epochs. The final parame¬ 
ters were estimated without using these SNe. The coefficients of 
the best-fit polynomials, R(y_x), ^(y-x) and the weighted stan¬ 
dard deviation around the fit are shown in Table [T^ They were 
derived from 22 LCs of 20 SNe with AM 15 and E(B - V) in the 
ranges 0.79+1.55 mag and -0.05 + +0.16 mag, respectively. The 
corrected color curves up to phases +45 days and the best fits are 
shown in Fig.[T^ for clarity omitting the error bars. 

One can see that the derived R(y-x) parameters are smaller 
than what can be expected. From the Hubble diagram analyses 
we estimated Ry ^ 1.9 and Rj/h - 0.3, and thus R(y-x) should 
be around 1.6. However, given the small range of E{B - V) and 
the relatively large scatter, the sensitivity to R(y-x) is not very 

'* this was done during the LC fitting. 


Article number, page 20 of 


38 


good. We also note that the primary goal was to derive the un¬ 
reddened loci, the scatter around them and the dependence on 
AMi 5 for the purpose of estimating the V - J and V - H color 
excesses between phases -11 and +12 days. Ideally, in deriving 
the loci one should use individual R(y-x) parameters for each 
SN, e.g. as in Burns et ar](2014|, but the use of common value 


for these low-reddened SNe should have small effect on the final 
result. 

The scatter around the loci is 0.10-0.12 mag. This is slightly 
larger that the intrinsic scatter of V - J and V - H estimated by 


Burns et al. (2014 1 . However, we note that part of the scatter of 


the y-NIR colors in our analysis is not statistical. A closer ex¬ 
amination of the data shows that most SNe are either bluer or 
redder than the mean loci at all epochs. This behavior is likely 
to be the result of several factors. Using an average R(y-x) value 
instead of individual ones, systematic photometry calibration er¬ 
rors and intrinsic color scatter can all cause systematic offset of 
the color curves. The fact that four SNe from the low-reddened 
sample had systematically too blue colors clearly demonstrate 
this point. We note however that the offsets of SNe 20061e and 
20061f are too big to be explained by a different extinction law. 


5. Discussion 

5. 1. SNe la in the optical 

We analyzed the B and V LCs of the largest nearby SN la compi¬ 
lation to date - 474 LCs from 11 different samples. From a sam¬ 
ple of 348 SNe with redshift above z = 0.015 we found that the 
weighted scatter around the best-fit Hubble line is ~ 0.18 + 0.01 
mag, independent on the assumed nature of the intrinsic scat¬ 
ter (color or luminosity). However, the scatter of the individual 
samples differs significantly (Table [T4)l and it is between 0.14 
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Table 16. Parameters of the 3rd-order polynomials describing the V — J 
and V — H colors evolution and the R(v-x) and a(v-x) parameters. 


Parameter 

V-J 

V-H 

Co 

-0.59378 

-0.83613 

Cl 

-0.07383 

-0.05885 

C2 

-2.66523E-03 

2.88438E-03 

C3 

-1.17256E-07 

1.94898E-04 

n(v~X) 

0.320 (0.040) 

0.518 (0.052) 

R(V-X) 

0.897 (0.126) 

0.800 (0.150) 

O" intr 

0.095 

0.106 

RMS [-11,-h12] 

0.100 

0.116 


and 0.30 mag for the four largest samples (CSP, Berkeley, CfA3 
and CfA4). The corresponding intrinsic scatter is also different. 
These four samples have similar redshift distributions and photo¬ 
metric quality, and the scatter that can be introduced in the Hub¬ 
ble diagram by the galaxies peculiar velocities and the photomet¬ 
ric uncertainties should be similar. The fact that we find differ¬ 
ent cTina indicates the presence of additional sample-dependent 
source of uncertainty beyond the possible intrinsic SN la lumi¬ 
nosity scatter because the latter should make the same contribu¬ 
tion in all samples. 

In the color scatter model the smallest intrinsic scatter ixintr = 
0.1 mag is found in the CSP and CfA3 samples. From a similarly 
heterogeneous low-redshift sample, Betoule et al.| ( |2014| find an 
intrinsic luminosity scatter ~ 0.12 mag, but without including 
color scatter in their analysis. In our luminosity scatter model 
only the CSP sample has a comparable intrinsic scatter ~ 0.13 
mag and in the other big samples it is larger ~ 0.18 mag. In 
the smaller samples is within the range found the four big 
samples; in some it even approaches the small value of the CSP 
sample. Our new sample fits well in the middle of the range. 

From Table 14 it is seen that the mean absolute magnitudes 
vary by up to 0.16 mag between the samples. These relatively 
large differences are somehow surprising because when compar¬ 
ing the peak magnitudes of the SNe in common in the largest SN 
samples we found differences only up to 0.05 mag (Tablej^. The 
absolute magnitudes are derived from the full samples, while for 
the peak magnitude comparison only the SNe in common are 
used. Nevertheless, Berkeley and CfA3 samples have 50 SNe in 
common, and CSP and CfA3 30. If the difference of the abso¬ 
lute magnitude between these samples reflect systematic uncer¬ 
tainties in the photometry, it is reasonable to expect that these 
uncertainties will also be noticeable when comparing the peak 
magnitudes of 30-50 SNe in common. Whatever the reason for 
the differences of the estimated absolute magnitudes is, possible 
calibration and other systematic uncertainties are of serious con¬ 
cern for cosmology and Betoule et al. (20141 have demonstrated 
how important is that these cross-sample mismatches are under 
control. 

Analyzing the distribution of the maximum B - V colors of 
393 SNe la we found evidence that the intrinsic B - V color 
of normal SNe la at maximum light is a non-linear function of 
AMi 5 . Using our (B - - AM 15 relation given by Eq.l^to 

estimate E(B - y)obs, we find that the nearby SN la data is ^st 
described with Rb parameter ~ 2.9 or ~ 2.7 depending on the 
assumed origin of the intrinsic scatter. [Scolnic et 3T] ( |2014| l have 
suggested that by treating the color scatter properly one should 
recover the typical Milky Way value of 7?^ = 4.1. Because of the 
large span of the colors of the nearby SN la sample our estima¬ 
tion of Rb ~ 2.9 is robust and it is easy to see from Fig. 11 that 
Rb -4.1 is inconsistent with the data. 


5.2. SNe la in the NIR 

Despite the great success of the optical observations of SNe la, 
the accuracy of the derived cosmological parameters is still lim¬ 
ited by a number of systematic uncertainties (see, e.g., Conley 
|et al.|| 20 TT| [Betoule et al.||20l4) l. Two sources of such uncer¬ 
tainties are the assumptions that the slope of the "luminosity - 
light curve shape" relation and the properties of the dust grains 
dimming the SN light do not change over cosmic time, i.e. the 
ax and Rx parameters in Eq. are the same for all SN^^ The 
assumption for constant Rx is likely not correct and there is in¬ 
creasing evidence for its variation between different events in the 
nearby universe (Folatelli et al. 2010[ Burns et al. 2014[ Amanul- 
jlah et al.|2015| and references therein). There are also indications 
that the most heavily reddened SNe tend to have low Ry ~ 1-2, 
while in the more moderately reddened events a range of values 
up to Ry ~ 3 is observed. In addition, we do know that the dust 
properties in the Milky Way vary depending on the line of sight 
and this is likely to be the case for other galaxies too. Even if 
the reddening law were to remain unchanged, it is already prob¬ 
lematic that different populations would be observed as the red- 
shift increases. Eigurel^shows that the distribution of (B- V)b„,„ 
color of the nearby SNUa sample has an extended red tail due to 
dust extinction. With increasing the redshift more and more of 
the reddened SNe will be missed and the high-redshift sample 
will be dominated by low-reddening events. Even if one is to ap¬ 
ply color cut to the nearby sample, the distribution of the colors 
will still be different between the low- and high-redshift sample, 
which can introduce systematic differences. 

These difficulties are clearly mitigated when shifting the ob¬ 
servations to the rest-frame NIR. We find Rjjh - 0.2 - 0.4 and 
nj/H - 0.15 - 0.20. Compared to the optical, the correction for 
the LC shape is by a factor of at least 3.5 smaller and the extinc¬ 
tion corrections are by a factor of ~ 5 - 15 smaller. Thus, the ef¬ 
fect of possible systematic uncertainties associated with the LC 
shape and the extinction corrections are greatly reduced in the 
NIR. Eurther, our analysis of the largest compilation of SN la 
NIR LCs to date confirms previous conclusions that SNe la are 
as good, if not better, standard candles in the NIR as in the op¬ 
tical. With corrections for the LC shape and extinction the ob¬ 
served scatter around the Hubble line becomes ~ 0.14 + 0.01 
mag, which is similar to the best published value of 0.14 mag 
derived by [Betoule et al.j \20\4\ from optical observations of 
nearby SNe la. After applying the LC shape and the extinction 
corrections it was noticed that the Hubble residuals showed a 
correlation with the Tmax - 7/max color in both J and H. Adding 
the 7max - //max term in the regression analysis leads to further 
small reduction of the scatter of the whole sample in J band, but 
a slight increase is observed in H (Table [T5]l. Even though our 
analysis indicated that this term is statistically significant, fur¬ 
ther observations are needed to confirm this finding. It should 
also be noted that while previous studies have claimed that the 
H band Hubble diagram has considerably smaller scatter than J 
(e.g., Wood-Vasey et al.|2008 1 , we find that the scatter is equal 
in J and H. 

More importantly, our analysis suggests that the NIR intrin¬ 
sic luminosity scatter of SNe la may be better than in the op¬ 
tical. Estimating this quantity is a difficult task because it re¬ 
quires identifying all other sources of uncertainty and correctly 
accounting for them. In addition, low-redshift samples are sen¬ 
sitive to the assumed value of the galaxy velocity dispersion. 


Small evolution of the mean SN la luminosity driven by some global 
parameter such as the mean metallicity of the universe cannot be ex¬ 
cluded also, but this would be very difficult to probe observationally. 
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the exact value of which is rather uncertain. Considering that 
this latter uncertainty affects greatly the estimation of crjnti in 
the more nearby samples and the quality of the data in the two 
high-redshift samples, it appears that the BN12 sample is best 
suited to estimate crintr in the NIR. Even without applying any 
corrections crjntr in BN12 sample is rather small, 0.14 mag in J 
and only 0.074 mag in H (Table [BJ. Including corrections for 
the LC shape and the extinction further reduced cTjntr and in the 
BN12 sample it becomes 0.11 mag in J and 0.062 mag in H. 
Finally, adding a correction term proportional to - T/max iy- 
term) in the regression analysis leads to further decrease of crintr 
to 0.028 mag in J band for BN12, but a slight increase in H . 
The y-term also affects the estimation of in J band in our 
sample, reducing it from 0.23 mag to 0.08 mag. However, con¬ 
sidering that the y-term does not have the same big effect on ixintr 
in J band of the other samples, the decrease that we see in our 
and BN 12 samples may be a statistical fluctuation cause by the 
somehow small sample sizes. 

The particularly small cri„tr = 0.028 mag in J band of BN12 
sample is interesting. However, we have traced this small value 
to the fact that we included estimation of the uncertainties of the 
/T-correction in the error budget. The combination of excellent 
LC quality, redshift range and size of the ^T-corrections uncer¬ 
tainty makes the BN 12 sample most sensitive to the /T-correction 
uncertainty. We have assumed that the scatter of the individual 
/T-corrections reflects real difference on the energy distribution 
in the NIR spectra of different SNe, hence this scatter directly 
adds to the uncertainty of the peak magnitudes. However, this 
may not necessarily be the case. It is possible that (at least part 
of) the observed scatter is introduced by our procedure to pre¬ 
pare the spectra in the form suitable for computing /T-corrections 
(Sec. |3.3| l. Matching the spectra to photometry, which potentially 
may contain calibration errors, the interpolation across the tel¬ 
luric band, etc., certainly introduce extra scatter in the individual 
/T-correction. However, when using smoothing-spline interpola¬ 
tion to compute the /T-corrections at the phases of the observa¬ 
tion, this additional scatter should be smoothed out and should 
not be added to the peak magnitudes. More high-quality photom¬ 
etry and spectroscopy of nearby SNe la are needed in order to 
separate the intrinsic variations from the other sources of scatter 
in the ^T-corrections. For these reasons in Tables [MTTI the uncer¬ 
tainties of the peak magnitudes estimated by the LC fitter and the 
uncertainties of the /T-corrections are given separately. Without 
addition the uncertainty of the /T-corrections crjntr = 0.06 mag 
in J band in BN12 sample and the values of o-i^a in the other 
samples/bands are only weakly affected. Given all the uncertain¬ 
ties, we can conclude that the NIR intrinsic luminosity scatter of 
SNe la is likely to be around 0.08-0.10 mag. This is smaller than 
both, our estimate from the optical observations and the value of 
0.12 mag derived by |Betoule et al.| ( |2014| l with the SALT2 LC 
fitter. 


5.3. NIR observations at high redshift 

In order to provide constrains of the cosmological parameters 
with rest-frame NIR observations of SNe la, a large sample of 
events out to redshift z = 1 needs to be obtained. However, this 
will be very demanding even for the future biggest telescopes. 
Therefore the main goal of this study was to test the accuracy 
of estimating the NIR peak magnitudes of SNe la with only a 
single observation coupled with good optical data, which is the 
most efficient approach to use rest-frame NIR observations of 
SNe la in cosmology. For our purposes we divided the published 
SNe in three groups GOLD, SILVER and BRONZE according to 
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the sampling of their NIR LCs and analyzed them together with 
our new sample of 16 SNe and the BN12 sample, both of which 
are at higher redshift. 

The analysis of the GOLD sample and a series of simulations 
gave strong indications that in principle it should be possible to 
estimate accurately the NIR peak magnitudes with a single ob¬ 
servations within 10 days from B band maximum light. By de¬ 
sign our new sample contains SNe with single NIR observations 
and only two SNe have two observations. Out of 28 SNe in the 
BRONZE sample, 16 have at most two observations per band, 
and all the others have three or more observations in at least one 
band. Out of the 15 SNe in the SILVER sample, only 4 have 
2 or less observations. Thus, along with our new sample, the 
BRONZE one is also well-suited to test the methodology. 

From the analysis of the two best samples, GOLD and BN12, 
we concluded that the intrinsic luminosity scatter ixintr of SNe la 
in the NIR is likely around 0.08-0.10 mag. We find that cTintr in¬ 
ferred from the other three samples is larger. Although, we have 
identified and took into account a number of uncertainties in¬ 
troduced when using under-sampled NIR LCs, it appears likely 
that some additional source(s) of scatter is also present. We es¬ 
timate that this additional scatter should be at the level of ~ 0.1 
mag in all cases, except for our sample in J with no y-term cor¬ 
rection applied, where it could be 0.2 mag. However, after in¬ 
cluding the y-term cri„tr in J band in our sample also becomes 
rather small 0.08 mag. In addition, we found that the absolute 
magnitudes estimated from the GOLD, SILVER, BRONZE and 
our new samples are in a good agreement, indicating that using 
under-sampled LCs does not lead to significant bias in the esti¬ 
mated NIR peak magnitudes. 

Well-sampled NIR LCs provide the best estimate of the peak 
magnitudes and the smallest scatter on the Hubble diagram. 
However, such LCs will be very difficult, if possible at all, to 
obtain at redshifts beyond z ~ 0.5 with current ground-based fa¬ 
cilities. NIR observations of SNe la beyond z ~ 0.5 will only 
be possible with the next generation of 40-m class telescopes 
(ESO Extremely Large Telescope, The Thirty Meter Telescope) 
and JWST In our discussion we focus on JWST because its ex¬ 
posure time calculator (ETC) is available and we can compute 
the exposure times needed. In the redshift range z = 0.6 -1.2 
the rest-frame J and H bands are shifted to the JWST F200W, 
F277W and F356W bands, respective!}!*^ The observed peak 
(Vega) magnitudes of a typical SN la at these bands will be 
~ 23-24 mag. The JWST ETC estimates that with 15 min expo¬ 
sure it should be possible to obtain photometry with S/N^ 20-80 
at NIR maximum light and with S/N^ 10 - 45 at 10 days past 
tB^. At maximum, observations out to z ~ 1.8 will also be pos¬ 
sible with exposure times of 40-60 min, but at z > 1.3 the rest- 
frame optical observations and the SN discovery also become 
problematic as it needs to be dome in the NIR from space-based 
facilities. 

However, obtaining LCs with five points (plus one more if 
galaxy reference image is needed) is a considerable investment 
in terms of observing time at the world’s most advanced and 
expensive observatory. Our results suggest that a better strategy 
may be to obtain a larger SN sample with single observation, 
preferably within 2-3 days from the NIR maximum. Depending 
on whether galaxy reference images are needed, with this ap¬ 
proach it will be possible to observe by a factor of 3 to 5 more 
SNe. Considering the values of the intrinsic luminosity scatter 
that we estimate in the different samples and the results of the 


note however that the JWST bands are twice broader than the red- 
shifted J and H. 
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simulations, this appears advantageous. The SN discoveries will 
likely come from LSST, which will also have the capability to 
provide B and V LCs out to redshift z ~ 0.9 ( Ivezic et al.|2008| l. 
Taking into account that the minimum technical time to trigger 
JWST will likely be around 2 days, to observe a SN within a few 
days from NIR maximum, the SN candidates need to be spec¬ 
troscopically classified as SNe la at least one week before 
Very effective selection of SNe la for spectroscopic screening 
(90% success) can be made by studying the early multi-color 
LCs of the candidates, as shown by |Sako et al. ( |2008| l. Thus, to¬ 
gether with a dedicated program to classify SNe la on 8 -m class 
telescope LSST should provide enough targets for NIR follow¬ 
up. 


6. Conclusions 


We obtained multi-epoch optical UBVRI and single-epoch NIR 
J and H photometry of 16 SNe la in the redshift range z - 
0.037 - 0.183. This sample doubles the leverage of the current 
SN la NIR Hubble diagram. We analyzed this new data set to¬ 
gether with 102 NIR and 458 optical LCs of normal SNe la from 
the literature. Our findings can be summarized as follows. 

• a sub-sample of SNe with well-sampled NIR LCs (GOLD 
sample) allows template fitting with the stretch and the time of 
maximum as free parameters. We found that the optical and NIR 
stretches are nearly one-to-one related and that the time of B- 
band maximum light estimated from the NIR and optical LC fits 
are very close. This further reinforces the results of |Krisciunas 
et al. ( |2004c|b| l and |Wood-Vasey et ah ( |2008| l who showed that 
smooth LC templates, the time axis of which is stretched accord¬ 
ing to the optical stretch parameter, accurately describe the rest- 
frame J and //-band light curves within 10 days from B-band 
maximum. This allowed us to estimate the NIR peak magnitudes 
of the SNe with under-sampled LCs, including our new sample. 

• Our analysis confirms previous findings that SNe la are 
at least as good standard candles in the NIR as in the optical, 
now with good statistics of > 100 objects. We find evidence that 
the NIR Hubble residuals weakly correlate with the LC shape 
parameter AM 15 and the color excess E(B - y)obs, and for the 
first time we report a possible dependence on the Tmax - //max 
color. The intrinsic luminosity scatter of SNe la in the NIR is 
estimated to be around 0.08-0.10 mag. The weighted rms scatter 
of the Hubble diagram is 0.14 + 0.01 mag in both J and H, and 
SNe la appear to be equally good standard candles in both bands. 

• compared to the optical, the corrections for the LC shape in 
the NIR are by a factor of 3.5 smaller. The extinction corrections 
are also smaller by an order of magnitude and consistent with 
low Ry - 1.9. Thus, possible systematic uncertainties associ¬ 
ated with the LC shape and the extinction corrections are greatly 
reduced in the NIR. 

• We find that the intrinsic B - V color of SNe la at B-band 

maximum light, (B - can be described by the following 

non-linear function of AM 15 : 


(B - = 0.036 exp 


AMi 5 - 1.3 


0.3455 


0.084(+0.01). 


The intrinsic scatter of (B - V)b„„,o was estimated to be ~ 0.04 
mag. 

• the weighted scatter of the B band Hubble residuals is 
0.184 + 0.007 mag, independent on the assumed nature of the 
intrinsic scatter (color or luminosity). The inferred mean abso¬ 
lute magnitudes and intrinsic luminosity scatter was found to 
vary between the samples, indicating possible systematic uncer¬ 
tainties. The dust extinction in the host galaxies is on average 


described with Ry parameter ~ 1.9 or ~ 1.7 depending on the 
assumed origin of the intrinsic color scatter. Because of the large 
span of the maximum B - V color of the nearby SN la sample 
our estimation of Ry is robust and Ry = 3.1 is inconsistent with 
the data. 

In conclusion, we have now extended the NIR SN la Hubble 
diagram to its non-linear part at z ~ 0 . 2 , i.e. where the lever 
arm in redshift makes it possible to start probing the cosmic 
constituents, most notably the properties of dark energy. Inter¬ 
estingly, our study confirms that it is feasible to accomplish this 
result with very modest sampling of the NIR LCs, if comple¬ 
mented by well-sampled optical LCs to determine the SN phase 
and the optical stretch parameter. This bodes well for the future, 
in particular when LSST and JWST will be operational. JWST 
will have the capabilities to observe the rest-frame J and H out to 
redshift z - 1.8 and LSST will provide well-sampled rest-frame 
B and V LCs out to redshift z - 0.9. In the view of our results, 
we suggest that the most efficient way to extend the SN la NIR 
Hubble diagram to high redshifts will be to observed a large of 
SNe la with a single observation close to the NIR maximum. 
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Table 4. S-corrected SN photometry. The number in parentheses are the uncertainties. 


JD-2400000 

U 

B 

y 

R 

/ 

Instrument 

SN 2008bz 

54586.409 

18.499 (0.032) 

18.284 (0.033) 

18.182 (0.029) 

18.169 (0.019) 

18.765 (0.030) 

ALFOSC 

54592.441 

19.425 (0.055) 

18.935 (0.015) 

18.499 (0.014) 

18.596 (0.029) 

19.188 (0.077) 

MOSCA 

54598.436 

20.284 (0.069) 

19.631 (0.019) 

18.930 (0.012) 

18.766 (0.013) 

19.308 (0.024) 

ALFOSC 

54610.476 

- 

20.803 (0.042) 

19.739 (0.036) 

19.177 (0.029) 

19.139 (0.032) 

STANCam 

SNF20080510-005 

54599.395 

18.583 (0.033) 

19.055 (0.022) 

18.996 (0.020) 

18.927 (0.015) 

19.166 (0.022) 

ALFOSC 

54610.463 

19.206 (0.097) 

19.139 (0.036) 

18.875 (0.047) 

18.781 (0.019) 

19.246 (0.040) 

STANCam 

54618.569 

20.201 (0.058) 

19.644 (0.025) 

19.219 (0.021) 

19.186 (0.020) 

19.940 (0.044) 

ALFOSC 

54624.408 

20.837 (0.121) 

20.155 (0.029) 

19.577 (0.022) 

19.421 (0.020) 

20.154 (0.052) 

ALFOSC 

54636.493 

- 

21.393 (0.075) 

20.307 (0.043) 

19.782 (0.058) 

19.745 (0.034) 

ALFOSC 

SNF20080512-008 

54609.616 

19.734 (0.057) 

19.315 (0.032) 

18.851 (0.018) 

18.747 (0.019) 

19.306 (0.029) 

STANCam 

54613.610 

- 

- 

19.061 (0.037) 

18.931 (0.013) 

- 

ALFOSC 

54618.612 

20.953 (0.070) 

20.175 (0.018) 

19.330 (0.013) 

19.070 (0.013) 

19.519(0.023) 

ALFOSC 

54624.443 

21.133 (0.146) 

20.699 (0.037) 

19.617(0.016) 

19.145 (0.013) 

19.369 (0.020) 

ALFOSC 

54635.675 

- 

21.033 (0.146) 

20.219 (0.067) 

19.629 (0.027) 

19.559 (0.046) 

STANCam 

54641.549 

- 

21.521 (0.127) 

20.515 (0.055) 

19.918 (0.037) 

19.905 (0.032) 

ALFOSC 

SNF20080512-010 

54612.562 

19.413 (0.057) 

18.966 (0.025) 

18.568 (0.013) 

18.545 (0.013) 

19.197 (0.029) 

ALFOSC 

54618.597 

20.396 (0.050) 

19.843 (0.018) 

19.033 (0.012) 

18.784 (0.013) 

19.116(0.020) 

ALFOSC 

54626.470 

21.239 (0.111) 

20.813 (0.028) 

19.751 (0.017) 

19.111 (0.015) 

18.978 (0.016) 

ALFOSC 

54635.658 

- 

21.047 (0.168) 

20.279 (0.070) 

19.772 (0.040) 

19.602 (0.048) 

STANCam 

54641.523 

- 

21.414(0.115) 

20.560 (0.053) 

20.000 (0.036) 

19.947 (0.050) 

ALFOSC 

SNF20080516-000 

54611.490 

- 

18.422 (0.047) 

18.261 (0.020) 

18.179 (0.013) 

18.595 (0.022) 

STANCam 

54618.507 

19.114 (0.035) 

18.813 (0.012) 

18.537 (0.012) 

18.473 (0.011) 

19.086 (0.023) 

ALFOSC 

54626.425 

20.322 (0.073) 

19.628 (0.014) 

19.019(0.011) 

18.827 (0.012) 

19.321 (0.022) 

ALFOSC 

54635.401 

- 

20.529 (0.104) 

19.535 (0.041) 

19.032 (0.023) 

19.216 (0.099) 

STANCam 

54641.401 

- 

20.897 (0.059) 

19.873 (0.024) 

19.173 (0.014) 

19.141 (0.026) 

ALFOSC 

54656.432 

- 

21.256 (0.100) 

20.523 (0.048) 

19.936 (0.023) 

20.021 (0.056) 

ALFOSC 

SNF20080516-022 

54611.509 

18.229 (0.033) 

18.469 (0.019) 

18.393 (0.024) 

18.321 (0.027) 

18.751 (0.027) 

STANCam 

54618.583 

19.278 (0.032) 

18.895 (0.011) 

18.620 (0.012) 

18.573 (0.014) 

19.127 (0.019) 

ALFOSC 

54624.459 

20.013 (0.062) 

19.498 (0.018) 

18.980 (0.013) 

18.889 (0.014) 

19.502 (0.025) 

ALFOSC 

54636.515 

- 

20.820 (0.072) 

19.780 (0.045) 

19.224 (0.019) 

19.257 (0.030) 

ALFOSC 

54641.497 

- 

21.166 (0.050) 

20.235 (0.032) 

19.358 (0.023) 

19.361 (0.029) 

ALFOSC 

54648.564 

- 

21.453 (0.147) 

- 

- 

- 

ALFOSC 

54656.455 

- 

21.691 (0.078) 

20.857 (0.039) 

20.170 (0.028) 

20.298 (0.068) 

ALFOSC 

SNF20080517-010 

54610.471 

- 

19.907 (0.046) 

19.603 (0.017) 

19.479 (0.029) 

19.615 (0.023) 

STANCam 

54618.532 

21.150 (0.101) 

20.601 (0.023) 

19.965 (0.020) 

19.810(0.019) 

20.239 (0.086) 

ALFOSC 

54626.443 

22.413 (0.258) 

21.496 (0.042) 

20.534 (0.022) 

20.132 (0.019) 

20.314(0.074) 

ALFOSC 

54636.451 

- 

- 

- 

- 

20.273 (0.068) 

ALFOSC 

54641.416 

- 

22.729 (0.152) 

21.460 (0.061) 

20.868 (0.034) 

20.460 (0.061) 

ALFOSC 

54648.460 

- 

22.946 (0.356) 

21.864 (0.130) 

21.261 (0.135) 

20.810(0.155) 

ALFOSC 

SNF20080522-000 

54618.556 

16.486 (0.021) 

17.147 (0.013) 

17.126 (0.012) 

17.134 (0.011) 

17.330 (0.012) 

ALFOSC 

54625.591 

- 

17.269 (0.011) 

17.130 (0.021) 

17.009 (0.033) 

17.409 (0.075) 

ANDICAM 

54626.457 

17.073 (0.024) 

17.267 (0.017) 

17.123 (0.009) 

17.029 (0.012) 

17.491 (0.014) 

ALFOSC 

54628.555 

- 

17.389 (0.022) 

17.161 (0.017) 

17.075 (0.018) 

17.495 (0.066) 

ANDICAM 

54631.544 

- 

17.520 (0.025) 

17.260 (0.016) 

17.241 (0.046) 

17.706 (0.138) 

ANDICAM 

54635.560 

- 

17.846 (0.033) 

17.462 (0.029) 

- 

- 

ANDICAM 

54636.478 

18.231 (0.050) 

17.904 (0.023) 

17.586 (0.028) 

17.602 (0.016) 

18.141 (0.022) 

ALFOSC 

54645.584 

- 

19.059 (0.037) 

18.106 (0.019) 

17.801 (0.022) 

17.903 (0.059) 

ANDICAM 

54650.536 

- 

19.538 (0.026) 

18.381 (0.016) 

17.916 (0.016) 

17.923 (0.035) 

ANDICAM 
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Table 4. continued. 


JD-2400000 

U 

B 

y 

R 

/ 

Instrument 

54656.444 

- 

19.818 (0.018) 

18.662 (0.014) 

18.079 (0.011) 

18.076 (0.012) 

ALFOSC 

54658.547 

- 

19.929 (0.103) 

18.770 (0.049) 

18.172 (0.051) 

18.138(0.108) 

ANDICAM 

54663.537 

- 

20.190 (0.205) 

19.018 (0.065) 

18.438 (0.044) 

18.378 (0.067) 

ANDICAM 

54664.409 

- 

- 

18.980 (0.047) 

18.487 (0.021) 

18.447 (0.051) 

STANCam 

54671.400 

- 

20.303 (0.019) 

19.305 (0.013) 

18.805 (0.014) 

18.854 (0.015) 

ALFOSC 

SNF20080522-001 

54618.492 

18.926 (0.046) 

18.460 (0.015) 

18.052 (0.011) 

17.984 (0.012) 

18.608 (0.028) 

ALFOSC 

54626.412 

19.892 (0.051) 

19.308 (0.017) 

18.419(0.011) 

18.211 (0.009) 

18.567 (0.018) 

ALFOSC 

54639.589 

- 

20.195 (0.161) 

19.183 (0.129) 

18.427 (0.066) 

18.429 (0.106) 

ANDICAM 

54640.606 

- 

20.230 (0.124) 

19.212 (0.062) 

18.643 (0.051) 

18.520 (0.070) 

ANDICAM 

54651.522 

- 

20.643 (0.044) 

19.781 (0.045) 

19.104(0.045) 

19.204(0.106) 

ANDICAM 

54656.397 

- 

20.702 (0.126) 

19.785 (0.067) 

19.270 (0.034) 

19.407 (0.063) 

ALFOSC 

54661.531 

- 

20.870 (0.224) 

20.049 (0.141) 

19.385 (0.105) 

19.715 (0.200) 

ANDICAM 

54671.388 

- 

21.076 (0.099) 

20.200 (0.037) 

19.820 (0.034) 

19.976 (0.071) 

ALFOSC 

SNF20080522-004 

54618.475 

20.326(0.131) 

19.740 (0.031) 

19.529 (0.023) 

19.383 (0.026) 

19.932 (0.096) 

ALFOSC 

54626.398 

- 

20.364 (0.030) 

19.921 (0.023) 

19.762 (0.026) 

20.294 (0.078) 

ALFOSC 

54636.464 

- 

- 

20.635 (0.076) 

20.098 (0.045) 

20.174 (0.058) 

ALFOSC 

54641.388 

- 

21.563 (0.141) 

20.820 (0.067) 

20.263 (0.045) 

20.292 (0.073) 

ALFOSC 

54648.396 

- 

22.009 (0.558) 

21.151 (0.195) 

20.317 (0.115) 

- 

ALFOSC 

54656.386 

- 

- 

- 

20.913 (0.161) 

20.637 (0.125) 

ALFOSC 

SNF20080522-011 

54613.600 

- 

16.845 (0.012) 

16.888 (0.012) 

16.835 (0.012) 

17.061 (0.024) 

ALFOSC 

54624.393 

17.182 (0.023) 

17.085 (0.015) 

16.954 (0.012) 

16.911 (0.013) 

17.431 (0.016) 

ALFOSC 

54638.723 

- 

18.458 (0.068) 

17.800 (0.074) 

17.501 (0.032) 

17.763 (0.044) 

ANDICAM 

54650.585 

- 

19.556 (0.043) 

18.386 (0.024) 

17.819 (0.027) 

17.702 (0.034) 

ANDICAM 

54656.468 

- 

19.720 (0.022) 

18.717(0.014) 

18.194(0.011) 

18.111 (0.024) 

ALFOSC 

54661.632 

- 

20.178 (0.209) 

18.971 (0.062) 

18.458 (0.075) 

18.460 (0.073) 

ANDICAM 

54664.419 

- 

19.942 (0.093) 

19.036 (0.057) 

18.571 (0.025) 

18.648 (0.105) 

STANCam 

54671.413 

- 

20.108 (0.019) 

19.208 (0.014) 

18.798 (0.012) 

18.850 (0.023) 

ALFOSC 

SNF20080606-012 

54636.530 

19.298 (0.064) 

18.944(0.023) 

18.824 (0.019) 

18.724 (0.019) 

18.954 (0.034) 

ALFOSC 

54641.511 

19.861 (0.067) 

19.337 (0.017) 

19.017(0.015) 

18.951 (0.014) 

19.283 (0.024) 

ALFOSC 

54648.501 

- 

20.154 (0.041) 

19.435 (0.025) 

- 

19.396 (0.075) 

ALFOSC 

54649.677 

- 

20.264 (0.064) 

19.519 (0.052) 

19.331 (0.060) 

19.669 (0.149) 

ANDICAM 

54654.664 

- 

20.965 (0.170) 

19.862 (0.065) 

19.415 (0.062) 

19.387 (0.141) 

ANDICAM 

54663.406 

- 

21.650 (0.110) 

20.532 (0.192) 

19.778 (0.089) 

19.558 (0.095) 

ALFOSC 

54671.427 

- 

21.901 (0.060) 

20.958 (0.036) 

20.337 (0.025) 

20.127 (0.033) 

ALFOSC 

SNF20080610-003 

54641.536 

20.012 (0.079) 

19.456 (0.021) 

19.129 (0.016) 

18.992 (0.015) 

19.548 (0.027) 

ALFOSC 

54648.441 

20.706 (0.167) 

20.013 (0.036) 

19.553 (0.026) 

19.321 (0.028) 

20.093 (0.138) 

ALFOSC 

54655.527 

- 

20.840 (0.022) 

19.998 (0.021) 

19.596 (0.016) 

19.903 (0.060) 

STANCam 

54663.391 

- 

- 

20.412(0.106) 

19.779 (0.071) 

19.883 (0.152) 

ALFOSC 

54671.443 

- 

21.962 (0.052) 

20.985 (0.031) 

20.217 (0.019) 

20.114(0.029) 

ALFOSC 
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Table 10. SILVER and BRONZE Supernova EC parameters 
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Appendix A: Second-order extinction 


The second-order extinction term in U and B bands arises from 
the fact that the transmission of Earth’s atmosphere is not con¬ 
stant but decreases rapidly at wavelengths below 5000A. The 
change of the transmission across the U and B bands is large 
enough to cause stars with different spectral energy distributions 
(SED) to experience slightly different extinction. Thus, the at¬ 
mospheric extinction in U and B bands is not constant but it is a 
function of the object SED. It can be parametrized as; 


ks = k'^ + k”{B-V) 
ku - k'fj + k'^ (U - B), 


(A.l) 

(A.2) 


where k'^g and k'ljg are the first-order and the second-order ex¬ 
tinction. I’he second-order extinction can be neglected in VRI. 
Unfortunately, in the majority of the photometric works on SNe 
la the second-order extinction term has been neglected. Here we 
discuss the uncertainties that are caused by this. 

The usual practice in SN photometry is to calibrate the mag¬ 
nitudes of the stars in the field of the SN in several photometric 
nights. This is usually done by deriving the parameters in Eqs.[^ 
by observing Landolt| ( 1992 2009| l standard fields that contain 
both blue and red stars in photometric conditions, and applying 
the equations to the field stars. Then in non-photometric nights 
the SN magnitudes are calibrated using the field stars as tetrary 
standards. Simple re-arrangement of Eqs.[2 taking B band as an 
example gives: 


B = b + (ctB - k”X)(B-V)-k'gX + zpB, (A.3) 


If kg is not fitted (or fixed to a reasonable value-0.03), one 

would not measure the true color term ctB but rather cts - kg X, 
i.e. the measured color-term will depend on the airmass of the 
observations. The second-order extinction in U and B has been 
measured and it is about -0.03. Using the spectrophotometric 
atlas of |Stritzinger et al.| ( |2005| l, various U and B filters, and 
atmospheric transmission measured at different observatories we 
estimated kg and kjj by means of synthetic photometry to be 

-0.03 as well. Therefore, two standard fields observed airmass 

X = 1 and 2 will give color terms that differ by ~ 0.03. 

Another problem arises because the B-V color index of SNe 
la changes from ~ -0.1 around maximum light to 1.3 about a 
month later. If a SN is observed at airmass 1.5 the correction 
for the second-order extinction at maximum and at -1-30 days 
will be ~ 0.0 and 0.058 mag, respectively. When the second- 
order extinction term is not included in the calibration, this will 
not be corrected for. While around maximum the effect will be 
negligible, 1-2 month after maximum the B magnitudes maybe 
as much as 0.06 mag in error. This will directly affect the B - 
V color index, which is often used to estimate the host galaxy 
extinction. Therefore, the estimation of the peak magnitudes can 
be affected by as much as 0.18 mag for Milky Way type dust 
with Ry = 3.1. 


Appendix B: S-corrections 

Appendix B. 1: Photometric systems responses 


To recover the effective system responses we follow the meth¬ 
ods described in [Stanishev et al.|( 2007|l. Th e relevant informa¬ 
tion was either taken from Stanishev et al.| ( |2007| l or from the 
corresponding observatories. Whenever possible measured filter 
transmission and CCD QEs were used; if unavailable, specifi¬ 
cations from the manufacturer were used. The total system re¬ 
sponses were computed by multiplying the filter transmissions 


by (a) the CCD quantum efficiency (QE), (b) the reflectivity of at 
least two aluminum surfaces, (a) a telluric absorption spectrum 
at airmass 1.3. The lens and window transmissions were not in¬ 
cluded because this information was unavailable. Note that the 
transmission of the Earth atmosphere was not included in the cal¬ 
culation of the effective U and B bands. The reason is that when 
computing the instrument color terms the second-order extinc¬ 
tion has been taken into account. The effect of the atmospheric 
transmission on the U and B bands is to shift the bands effective 
wavelength redward with increasing the airmass of the observa¬ 
tions. Taking into account the second-order extinction meant that 
the computed color terms correspond to observations above the 
atmosphere. 

Synthetic magnitudes were calculated from |Stritzinger et al.| 
( |2005) l spectrophotometry of Landolt standard stars, 

=-2.51ogU /f °‘(/l)/?“'(T)c/Tj . (B.l) 

The difference between the synthetic and the observed photome¬ 
try was fitted as a function of the observed color indices to com¬ 
pute synthetic color-terms e.g. for B we have; 

^std _ ^syn ^ _ ystd) ^ (g 2) 


Eor the y/?/-bands, the were close to the observed ones 

In some cases small differences exceeding the uncertainty 
were accounted for by shifting the filter transmissions until 
matched zp^"" is the zero point for synthetic photometry 
with the band 

As in Stanishev et al. ( |2007| l we found that the synthetic U 
and B bands were always too blue. This is, to some extent, to be 
expected because the neglected optical elements like lenses or 
windows, anti-reflection and other coatings will tend to reduce 
the system sensitivity shortward of ~4000 A. To account for the 
net effect of these uncertainties we modified the U and B bands 
by multiplying them with a Sigmoid function; 


F(A;Ao,A)^ 


1 


1 + exp(-(A - Ao)/A) 


(B.3) 


The parameters Aq and A were selected so that matched 
for both 17 and B. 

Eollowing |Stanishev et al.| ( |2007| l we also re-derived the stan¬ 
dard Johnson-Cousins system responses with the revised filter 
functions of Bessell & Murphy ( 2012[ ). This procedure also gave 
the zero points zp^" for synthetic photometry with the standard 
system responses R^“^. 


Appendix B.2: Computing the S-corrections 

The procedure that we used to calculate the S -corrections was 
the same as in [Stanishev et al.| ( |2007| l, but slightly modified to 
account for the fact that the U and B bands did not include the 
atmospheric transmission. The calibration correction from natu¬ 
ral magnitudes at airmass X, to standard magnitudes is 
computed by means of synthetic photometry; 

= (B.4) 

where S is the S -correction; 

5 = -2.5 log (f ff°\A)R^“^(A)dA) + zp'^" 

+2.5 log (/ ff°\A)R‘'^\A)p(Af dA) - zpT. (B.5) 
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Here is the photon spectral energy distributions (SED) 

of the SN, and are the responses of the natural and 

the standard systems, respectively, and p(A) is the transmission 
of the atmosphere at airmass X = 1 above the site of o bserv ation. 
Here zp^“ is the zero point obtained by fitting Eq. B.2 to the 
synthetic magnitudes mj" computed with the natural passband 
multiplied by the Earth’s atmosphere transmission at airmass X, 
e.g.: 


f 


-2.5 log I I ff°\A)R--\A)p{AfdA 


(B.6) 


Comparing Eq.|A.3|and Eq.|B.4|we see that the term ctg{B-V)- 


k'g X{B - y) is just an approximation of the S -correction: 


S ^ cteiB -V)- k” X(B - V). 


(B.7) 


By the way the zero points in Eq. |B.5| are derived, S ^ 0.0 for 
AOV stars with all color indexes zero, e.g. Vega. This ensures 
that for normal stars the synthetic 5-correction will work the 
same way as the linear color-term corrections. 

The best way to compute the S -corrections for a SN is to 
use spectrophotometry of the same object. If that is not possible 
the mean spectrum of SNe la should be used. In our work we 
have used Hsiao et al. ( |2007| l spectral template, corrected for the 
Milky Way reddening and the redshift of each SN. 
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SN2008bz 


SNF20080510-005 


^SNF20080512-b08. 


SNF20080512-010 


51^2008051 $-022 


^SNF20080516-000- 


SNF20080517-010 


SNF20JD80522-000 


SNF20080fe06-012 ' 


SNF20080522-Q01 


SNF20080522-004 


^SNF2Qp80522Tbl1 


Fig. 1. Finding charts of the SNe observed in this paper. The SNs are located at the intersection of the two ticks. The numbers mark the local 
sequence stars. In most cases the number is right-up of the star and in a few cases only it is left-down. 
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Table 3. Calibrated magnitudes of the local stars around the SNe. The number in parentheses are the uncertainties. 


No. 

U 

-^obs 

B 

^obs 

V 

^obs 

R 

^obs 

/ 

^obs 

SN 2008bz 

1 

15.755 (0.013) 

1 

15.814(0.014) 

1 

15.220 (0.011) 

1 

14.895 (0.012) 

1 

14.552 (0.012) 

1 

2 

16.029 (0.013) 

1 

16.101 (0.013) 

1 

15.587 (0.011) 

1 

15.278 (0.012) 

1 

14.973 (0.012) 

1 

3 

17.657 (0.018) 

1 

16.858 (0.016) 

1 

15.935 (0.012) 

1 

15.393 (0.013) 

1 

14.933 (0.012) 

1 

4 

19.522 (0.036) 

1 

18.324(0.020) 

1 

17.159 (0.012) 

1 

16.398 (0.014) 

1 

15.766 (0.013) 

1 

5 

20.519(0.077) 

1 

19.187 (0.020) 

1 

17.945 (0.012) 

1 

17.143 (0.013) 

1 

16.433 (0.015) 

1 

6 

19.023 (0.027) 

1 

18.719 (0.020) 

1 

17.973 (0.017) 

1 

17.531 (0.018) 

1 

17.115 (0.015) 

1 

7 

20.961 (0.108) 

1 

19.816 (0.033) 

1 

18.233 (0.015) 

1 

17.169 (0.016) 

1 

15.672 (0.016) 

1 

8 

18.942 (0.027) 

1 

19.028 (0.025) 

1 

18.479 (0.012) 

1 

18.095 (0.024) 

1 

17.750 (0.013) 

1 
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Table 3. continued. 


No. 

U 

^obs 

B 

^obs 

V 

^obs 

R 

^obs 

/ 

^obs 

9 

19.111 (0.030) 

1 

19.300 (0.013) 

1 

18.733 (0.017) 

1 

18.408 (0.029) 

1 

18.047 (0.039) 

1 

10 

19.065 (0.028) 

1 

19.364 (0.039) 

1 

18.972 (0.019) 

1 

18.701 (0.017) 

1 

18.415 (0.017) 

1 

11 

21.601 (0.191) 

1 

20.547 (0.023) 

1 

19.112(0.014) 

1 

18.176 (0.015) 

1 

16.968 (0.014) 

1 

12 

21.224(0.150) 

1 

20.482 (0.033) 

1 

19.175 (0.014) 

1 

18.349 (0.018) 

1 

17.533 (0.033) 

1 

13 

20.605 (0.095) 

1 

20.133 (0.021) 

1 

19.228 (0.015) 

1 

18.692 (0.014) 

1 

18.177 (0.031) 

1 

14 

- 

- 

20.778 (0.053) 

1 

19.382 (0.021) 

1 

18.341 (0.016) 

1 

16.808 (0.016) 

1 

15 

19.956 (0.058) 

1 

20.164 (0.012) 

1 

19.695 (0.049) 

1 

19.335 (0.026) 

1 

18.963 (0.037) 

1 

16 

20.148 (0.058) 

1 

20.305 (0.015) 

1 

19.737 (0.027) 

1 

19.366 (0.015) 

1 

19.071 (0.020) 

1 

SNF20080510-005 

1 

18.666 (0.023) 

3 

17.572 (0.028) 

3 

16.046 (0.010) 

3 

15.065 (0.015) 

2 

- 

- 

2 

17.311 (0.015) 

3 

17.316(0.031) 

3 

16.691 (0.015) 

3 

16.325 (0.022) 

3 

15.957 (0.023) 

3 

3 

17.619 (0.021) 

3 

17.931 (0.010) 

3 

17.450 (0.013) 

3 

17.148 (0.015) 

3 

16.786 (0.015) 

3 

4 

19.090 (0.100) 

3 

18.482 (0.022) 

3 

17.591 (0.017) 

3 

17.030 (0.018) 

3 

16.530 (0.014) 

3 

5 

17.808 (0.038) 

3 

18.155 (0.017) 

3 

17.693 (0.014) 

3 

17.390 (0.011) 

3 

17.055 (0.013) 

3 

6 

18.845 (0.078) 

3 

18.945 (0.028) 

3 

18.371 (0.016) 

3 

18.045 (0.013) 

3 

17.688 (0.022) 

3 

7 

- 

- 

20.019 (0.016) 

3 

18.674 (0.018) 

3 

17.780 (0.019) 

3 

16.989 (0.016) 

3 

8 

19.740 (0.031) 

3 

19.542 (0.045) 

3 

18.810(0.003) 

3 

18.352 (0.009) 

3 

17.910(0.019) 

3 

9 

19.483 (0.034) 

3 

19.623 (0.011) 

3 

19.028 (0.022) 

3 

18.689 (0.009) 

3 

18.385 (0.010) 

3 

10 

20.913 (0.196) 

3 

20.417 (0.030) 

3 

19.033 (0.012) 

3 

18.093 (0.007) 

3 

16.999 (0.020) 

3 

11 

- 

- 

20.831 (0.050) 

3 

19.140 (0.012) 

3 

18.049 (0.028) 

3 

16.410 (0.014) 

3 

12 

18.823 (0.046) 

3 

19.460 (0.022) 

3 

19.241 (0.038) 

3 

18.917 (0.022) 

3 

18.346 (0.019) 

3 

13 

- 

- 

21.079 (0.065) 

3 

19.630 (0.045) 

3 

18.634 (0.018) 

3 

17.238 (0.015) 

3 

14 

- 

- 

21.175 (0.048) 

3 

19.692 (0.018) 

3 

18.725 (0.021) 

3 

17.438 (0.013) 

3 

15 

20.126 (0.062) 

3 

20.510(0.023) 

3 

20.048 (0.002) 

3 

19.723 (0.039) 

3 

19.461 (0.027) 

3 

SNF20080512-008 

1 

15.949 (0.015) 

1 

15.909 (0.008) 

2 

15.208 (0.004) 

2 

14.836 (0.009) 

2 

- 

- 

2 

16.568 (0.017) 

1 

16.042 (0.011) 

2 

15.211 (0.005) 

2 

14.750 (0.012) 

1 

- 

- 

3 

16.984 (0.018) 

1 

16.403 (0.009) 

2 

15.528 (0.003) 

2 

15.026 (0.002) 

2 

- 

- 

4 

16.215 (0.015) 

1 

16.223 (0.002) 

2 

15.608 (0.007) 

2 

15.255 (0.006) 

2 

- 

- 

5 

16.944 (0.006) 

2 

16.608 (0.005) 

3 

15.845 (0.005) 

3 

15.428 (0.015) 

3 

15.073 (0.069) 

3 

6 

16.539 (0.016) 

1 

16.570 (0.033) 

2 

15.940 (0.017) 

2 

15.555 (0.015) 

2 

15.160 (0.011) 

1 

7 

17.317(0.015) 

2 

16.854 (0.005) 

3 

16.013 (0.009) 

3 

15.518 (0.010) 

3 

15.087 (0.012) 

1 

8 

16.180 (0.003) 

2 

16.457 (0.002) 

3 

16.069 (0.002) 

3 

15.799 (0.008) 

3 

15.507 (0.018) 

3 

9 

17.265 (0.018) 

2 

16.982 (0.010) 

3 

16.215 (0.007) 

3 

15.789 (0.003) 

3 

15.370 (0.007) 

2 

10 

17.046 (0.005) 

2 

16.953 (0.010) 

3 

16.256 (0.003) 

3 

15.878 (0.009) 

3 

15.505 (0.007) 

2 

11 

16.974 (0.013) 

2 

17.143 (0.008) 

3 

16.620 (0.010) 

3 

16.303 (0.000) 

3 

15.957 (0.008) 

3 

12 

19.299 (0.055) 

2 

18.087 (0.004) 

3 

16.730 (0.006) 

3 

15.856 (0.004) 

3 

15.052 (0.012) 

1 

13 

17.595 (0.009) 

2 

17.617 (0.009) 

3 

17.005 (0.009) 

3 

16.640 (0.009) 

3 

16.270 (0.006) 

3 

14 

17.638 (0.019) 

2 

17.628 (0.003) 

3 

17.017 (0.006) 

3 

16.653 (0.019) 

3 

16.298 (0.010) 

3 

15 

17.311 (0.002) 

2 

17.598 (0.004) 

3 

17.134 (0.002) 

3 

16.820 (0.007) 

3 

16.489 (0.003) 

3 

16 

17.988 (0.004) 

2 

18.057 (0.014) 

3 

17.441 (0.008) 

3 

17.081 (0.013) 

3 

16.699 (0.003) 

3 

17 

17.932 (0.006) 

2 

18.053 (0.006) 

3 

17.500 (0.001) 

3 

17.166 (0.008) 

3 

16.803 (0.005) 

3 

18 

17.931 (0.002) 

2 

18.087 (0.007) 

3 

17.560 (0.005) 

3 

17.247 (0.005) 

3 

16.918 (0.004) 

3 

19 

18.284 (0.018) 

2 

18.276 (0.007) 

3 

17.631 (0.005) 

3 

17.252 (0.008) 

3 

16.875 (0.006) 

3 

20 

20.284 (0.064) 

2 

19.110(0.010) 

3 

17.690 (0.010) 

3 

16.763 (0.014) 

3 

15.833 (0.001) 

3 

21 

19.759 (0.048) 

2 

18.834 (0.014) 

3 

17.797 (0.002) 

3 

17.173 (0.010) 

3 

16.625 (0.006) 

3 

22 

19.427 (0.009) 

2 

18.724(0.003) 

3 

17.808 (0.007) 

3 

17.262 (0.004) 

3 

16.785 (0.002) 

3 

23 

18.823 (0.051) 

2 

18.631 (0.006) 

3 

17.885 (0.009) 

3 

17.452 (0.005) 

3 

17.019 (0.005) 

3 

24 

18.259 (0.012) 

2 

18.400 (0.014) 

3 

17.949 (0.012) 

3 

17.634 (0.015) 

3 

17.309 (0.008) 

3 

25 

19.820 (0.012) 

2 

19.054 (0.011) 

3 

18.108 (0.008) 

3 

17.515 (0.010) 

3 

17.006 (0.019) 

3 

26 

18.823 (0.021) 

2 

18.799 (0.007) 

3 

18.141 (0.006) 

3 

17.779 (0.005) 

3 

17.397 (0.018) 

3 

27 

19.193 (0.028) 

2 

18.940 (0.005) 

3 

18.166 (0.006) 

3 

17.724 (0.012) 

3 

17.314(0.008) 

3 

28 

20.670 (0.117) 

2 

19.735 (0.023) 

3 

18.376 (0.004) 

3 

17.474 (0.011) 

3 

16.480 (0.002) 

3 

29 

19.002 (0.018) 

2 

19.108 (0.023) 

3 

18.498 (0.007) 

3 

18.144(0.006) 

3 

17.796 (0.010) 

3 

30 

18.764 (0.005) 

2 

18.987 (0.009) 

3 

18.519(0.011) 

3 

18.212(0.005) 

3 

17.873 (0.011) 

3 

31 

20.381 (0.040) 

2 

19.609 (0.010) 

3 

18.618 (0.009) 

3 

17.973 (0.013) 

3 

17.373 (0.010) 

3 

32 

19.282 (0.035) 

2 

19.322 (0.017) 

3 

18.672 (0.006) 

3 

18.292 (0.016) 

3 

17.907 (0.008) 

3 

33 

20.689 (0.075) 

2 

19.799 (0.022) 

3 

18.736 (0.006) 

3 

18.126 (0.009) 

3 

17.593 (0.008) 

3 

34 

19.532 (0.182) 

2 

19.480 (0.011) 

3 

18.794 (0.016) 

3 

18.396 (0.006) 

3 

17.990 (0.009) 

3 

35 

- 

- 

20.321 (0.005) 

3 

18.812(0.006) 

3 

17.835 (0.019) 

3 

16.710(0.005) 

3 
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Table 3. continued. 


No. 

U 

^obs 

B 

^obs 

V 

^obs 

R 

^obs 


I 

^obs 

36 


- 

- 

20.064 

(0.039) 

3 

18.836 

(0.007) 

3 

18.017 

(0.015) 

3 

17.295 

(0.007) 

3 

37 

20.194 

(0.039) 

2 

19.757 

(0.012) 

3 

18.916 

(0.017) 

3 

18.419 

(0.013) 

3 

17.933 

(0.003) 

3 

38 

19.465 

(0.087) 

2 

19.587 

(0.008) 

3 

18.960 

(0.020) 

3 

18.577 

(0.017) 

3 

18.199 

(0.014) 

3 

39 

19.184 

(0.015) 

2 

19.455 

(0.020) 

3 

18.975 

(0.001) 

3 

18.649 

(0.018) 

3 

18.332 

(0.015) 

3 

40 

20.594 

(0.035) 

2 

19.929 

(0.007) 

3 

18.993 

(0.016) 

3 

18.436 

(0.021) 

3 

17.920 

(0.020) 

3 

41 

20.297 

(0.049) 

2 

19.840 

(0.015) 

3 

19.004 

(0.006) 

3 

18.518 

(0.004) 

3 

18.038 

(0.013) 

3 

42 

20.044 

(0.091) 

2 

20.048 

(0.022) 

3 

19.340 

(0.004) 

3 

18.971 

(0.025) 

3 

18.573 

(0.015) 

3 

43 


- 

- 

21.181 

(0.122) 

3 

19.634 

(0.028) 

3 

18.602 

(0.005) 

3 

17.321 

(0.008) 

3 

44 

19.703 

(0.113) 

2 

20.046 

(0.014) 

3 

19.663 

(0.013) 

3 

19.399 

(0.012) 

3 

19.039 

(0.003) 

3 

45 

19.792 

(0.138) 

2 

20.136 

(0.020) 

3 

19.667 

(0.029) 

3 

19.380 

(0.005) 

3 

19.032 

(0.014) 

3 

46 

20.997 

(0.110) 

2 

20.570 

(0.053) 

3 

19.683 

(0.014) 

3 

19.122 

(0.013) 

3 

18.644 

(0.027) 

3 

47 

20.842 

(0.004) 

2 

20.635 

(0.039) 

3 

19.724 

(0.004) 

3 

19.211 

(0.027) 

3 

18.738 

(0.017) 

3 

48 


- 

- 

21.342 

(0.013) 

3 

19.797 

(0.032) 

3 

18.782 

(0.006) 

3 

17.300 

(0.006) 

3 

49 


- 

- 

21.013 

(0.051) 

3 

19.813 

(0.002) 

3 

19.030 

(0.017) 

3 

18.258 

(0.012) 

3 

50 

22.579 

(0.005) 

2 

21.085 

(0.115) 

3 

19.830 

(0.020) 

3 

19.030 

(0.016) 

3 

18.344 

(0.015) 

3 

51 


- 

- 

21.715 

(0.085) 

3 

20.197 

(0.020) 

3 

19.215 

(0.006) 

3 

17.927 

(0.011) 

3 

52 


- 

- 

21.787 

(0.122) 

3 

20.384 

(0.022) 

3 

19.367 

(0.019) 

3 

17.881 

(0.008) 

3 

53 


- 

- 

22.473 

(0.089) 

3 

20.747 

(0.048) 

3 

19.526 

(0.016) 

3 

17.655 

(0.016) 

3 

SNF20080512-010 

1 

16.985 

(0.024) 

2 

16.839 

(0.009) 

3 

16.138 

(0.009) 

3 

15.756 

(0.002) 

2 

15.361 

(0.013) 

1 

2 

19.479 

(0.005) 

2 

18.426 

(0.006) 

3 

16.933 

(0.006) 

3 

15.924 

(0.005) 

2 


- 

- 

3 

17.406 

(0.003) 

2 

17.607 

(0.009) 

3 

17.056 

(0.008) 

3 

16.733 

(0.009) 

3 

16.367 

(0.011) 

3 

4 

19.239 

(0.032) 

2 

18.354 

(0.010) 

3 

17.351 

(0.006) 

3 

16.738 

(0.005) 

3 

16.196 

(0.011) 

3 

5 

18.066 

(0.018) 

2 

18.071 

(0.036) 

3 

17.408 

(0.020) 

3 

17.037 

(0.012) 

3 

16.654 

(0.008) 

3 

6 

20.373 

(0.115) 

2 

19.086 

(0.028) 

3 

17.660 

(0.006) 

3 

16.722 

(0.011) 

3 

15.795 

(0.002) 

2 

7 

20.927 

(0.078) 

2 

19.750 

(0.010) 

3 

18.286 

(0.008) 

3 

17.329 

(0.014) 

3 

16.241 

(0.001) 

3 

8 

20.458 

(0.112) 

2 

19.528 

(0.004) 

3 

18.560 

(0.008) 

3 

17.934 

(0.001) 

3 

17.413 

(0.012) 

3 

9 

18.794 

(0.013) 

2 

19.087 

(0.009) 

3 

18.661 

(0.014) 

3 

18.365 

(0.013) 

3 

18.040 

(0.002) 

3 

10 

19.868 

(0.016) 

2 

19.688 

(0.004) 

3 

18.902 

(0.013) 

3 

18.441 

(0.008) 

3 

17.944 

(0.015) 

3 

11 

19.051 

(0.022) 

2 

19.407 

(0.015) 

3 

18.922 

(0.009) 

3 

18.612 

(0.013) 

3 

18.263 

(0.014) 

3 

12 

19.970 

(0.057) 

2 

19.881 

(0.035) 

3 

19.215 

(0.014) 

3 

18.829 

(0.010) 

3 

18.437 

(0.025) 

3 

13 

21.373 

(0.147) 

2 

20.356 

(0.009) 

3 

19.218 

(0.013) 

3 

18.518 

(0.010) 

3 

17.920 

(0.014) 

3 

14 


- 

- 

20.519 

(0.025) 

3 

19.261 

(0.019) 

3 

18.430 

(0.011) 

3 

17.739 

(0.012) 

3 

15 

19.461 

(0.021) 

2 

19.710 

(0.020) 

3 

19.282 

(0.008) 

3 

19.005 

(0.011) 

3 

18.695 

(0.027) 

3 

16 

19.854 

(0.048) 

2 

19.989 

(0.003) 

3 

19.418 

(0.017) 

3 

19.071 

(0.004) 

3 

18.697 

(0.031) 

3 

17 

19.567 

(0.008) 

2 

19.969 

(0.054) 

3 

19.460 

(0.011) 

3 

19.184 

(0.022) 

3 

18.848 

(0.024) 

3 

18 

20.183 

(0.148) 

2 

20.269 

(0.015) 

3 

19.516 

(0.009) 

3 

19.119 

(0.006) 

3 

18.712 

(0.009) 

3 

19 

19.687 

(0.019) 

2 

20.030 

(0.008) 

3 

19.562 

(0.007) 

3 

19.253 

(0.008) 

3 

18.910 

(0.048) 

3 

20 

22.052 

(0.158) 

2 

21.141 

(0.016) 

3 

19.622 

(0.037) 

3 

18.520 

(0.010) 

3 

16.812 

(0.006) 

3 

21 

22.098 

(0.015) 

2 

21.061 

(0.011) 

3 

19.687 

(0.006) 

3 

18.800 

(0.005) 

3 

17.942 

(0.010) 

3 

22 

20.662 

(0.038) 

2 

20.535 

(0.028) 

3 

19.772 

(0.002) 

3 

19.336 

(0.008) 

3 

18.918 

(0.024) 

3 

23 

21.740 

(0.140) 

2 

21.052 

(0.047) 

3 

19.891 

(0.028) 

3 

19.087 

(0.017) 

3 

18.378 

(0.009) 

3 

24 


- 

- 

21.424 

(0.078) 

3 

20.036 

(0.024) 

3 

19.110 

(0.012) 

3 

18.065 

(0.011) 

3 

25 


- 

- 

21.635 

(0.147) 

3 

20.276 

(0.026) 

3 

19.370 

(0.008) 

3 

18.407 

(0.015) 

3 

26 


- 

- 

21.997 

(0.067) 

3 

20.508 

(0.032) 

3 

19.451 

(0.010) 

3 

17.912 

(0.013) 

3 

27 

20.278 

(0.046) 

2 

20.853 

(0.042) 

3 

20.514 

(0.024) 

3 

20.111 

(0.028) 

3 

19.695 

(0.028) 

3 

SNF20080516-000 

1 


- 

- 

15.750 

(0.004) 

2 

15.145 

(0.003) 

2 


- 

- 


- 

- 

2 


- 

- 

16.011 

(0.000) 

2 

15.221 

(0.007) 

2 


- 

- 


- 

- 

3 

16.989 

(0.026) 

1 

16.337 

(0.011) 

3 

15.405 

(0.007) 

3 

14.899 

(0.016) 

1 


- 

- 

4 

16.639 

(0.022) 

1 

16.864 

(0.004) 

4 

16.312 

(0.006) 

4 

15.983 

(0.008) 

4 

15.629 

(0.013) 

3 

5 

18.985 

(0.046) 

1 

18.021 

(0.004) 

4 

16.876 

(0.008) 

4 

16.160 

(0.008) 

4 

15.562 

(0.012) 

3 

6 

17.486 

(0.024) 

1 

17.495 

(0.017) 

3 

16.909 

(0.007) 

4 

16.550 

(0.004) 

4 

16.202 

(0.012) 

4 

7 

18.227 

(0.031) 

1 

17.782 

(0.004) 

4 

16.914 

(0.003) 

4 

16.408 

(0.008) 

4 

15.951 

(0.015) 

4 

8 

19.533 

(0.067) 

1 

18.844 

(0.019) 

4 

17.777 

(0.004) 

4 

17.112 

(0.009) 

4 

16.547 

(0.011) 

4 

9 

20.700 

(0.153) 

1 

19.701 

(0.012) 

4 

18.069 

(0.014) 

4 

17.001 

(0.005) 

4 

15.506 

(0.021) 

4 

10 

20.644 

(0.146) 

1 

19.436 

(0.019) 

4 

18.106 

(0.007) 

4 

17.223 

(0.011) 

4 

16.288 

(0.014) 

4 

11 

18.760 

(0.041) 

1 

18.802 

(0.035) 

4 

18.172 

(0.009) 

4 

17.828 

(0.012) 

4 

17.442 

(0.022) 

4 

12 


- 

- 

19.641 

(0.025) 

4 

18.325 

(0.017) 

4 

17.482 

(0.010) 

4 

16.612 

(0.017) 

4 

13 

20.302 

(0.118) 

1 

19.659 

(0.020) 

4 

18.346 

(0.008) 

4 

17.493 

(0.017) 

4 

16.729 

(0.017) 

4 
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Table 3. continued. 


No. 

U 

^obs 

B 

^obs 

V 

^obs 

R 

^obs 


I 

^obs 

14 

18.990 

(0.044) 

1 

19.031 

(0.013) 

4 

18.347 

(0.010) 

4 

17.952 

(0.017) 

4 

17.549 

(0.017) 

4 

15 


- 

- 

21.106 

(0.050) 

4 

19.621 

(0.013) 

4 

18.539 

(0.020) 

4 

16.936 

(0.011) 

4 

SNF20080516-022 

1 

16.508 

(0.006) 

2 

16.486 

(0.005) 

4 

15.844 

(0.005) 

4 

15.468 

(0.008) 

3 


- 

- 

2 

16.712 

(0.002) 

2 

16.591 

(0.009) 

4 

15.884 

(0.003) 

4 

15.462 

(0.007) 

3 


- 

- 

3 

17.385 

(0.001) 

2 

17.370 

(0.001) 

4 

16.780 

(0.007) 

4 

16.425 

(0.007) 

4 

16.081 

(0.004) 

4 

4 

19.424 

(0.021) 

2 

18.286 

(0.003) 

4 

17.122 

(0.010) 

4 

16.388 

(0.011) 

4 

15.775 

(0.005) 

4 

5 

18.680 

(0.026) 

2 

18.456 

(0.003) 

4 

17.706 

(0.008) 

4 

17.274 

(0.017) 

4 

16.884 

(0.006) 

4 

6 

18.836 

(0.057) 

2 

19.086 

(0.013) 

4 

18.582 

(0.007) 

4 

18.274 

(0.011) 

4 

17.947 

(0.010) 

4 

7 

20.297 

(0.067) 

2 

19.775 

(0.008) 

4 

18.914 

(0.009) 

4 

18.367 

(0.014) 

4 

17.883 

(0.017) 

4 

8 

19.459 

(0.009) 

2 

19.758 

(0.027) 

4 

19.328 

(0.030) 

4 

19.045 

(0.013) 

4 

18.781 

(0.011) 

4 

9 

21.549 

(0.089) 

2 

20.678 

(0.026) 

4 

19.608 

(0.021) 

4 

18.929 

(0.021) 

4 

18.339 

(0.019) 

4 

10 


- 

- 

21.127 

(0.023) 

4 

19.656 

(0.017) 

4 

18.709 

(0.016) 

4 

17.479 

(0.011) 

4 

11 


- 

- 

21.659 

(0.076) 

4 

20.318 

(0.025) 

4 

19.373 

(0.018) 

4 

18.103 

(0.010) 

4 

SNF20080517-010 

1 

16.357 

(0.062) 

3 

16.562 

(0.009) 

4 

16.058 

(0.021) 

4 

15.744 

(0.009) 

4 

15.420 

(0.022) 

3 

2 

19.073 

(0.008) 

2 

19.381 

(0.023) 

3 

18.957 

(0.025) 

3 

18.677 

(0.007) 

3 

18.375 

(0.026) 

3 

3 

19.256 

(0.038) 

3 

20.335 

(0.025) 

4 

20.131 

(0.012) 

4 

19.755 

(0.022) 

4 

19.376 

(0.017) 

3 

4 


- 

- 

21.247 

(0.035) 

3 

20.150 

(0.025) 

3 

19.466 

(0.005) 

3 

18.874 

(0.014) 

3 

5 


- 

- 

21.679 

(0.114) 

4 

20.297 

(0.026) 

4 

19.251 

(0.027) 

4 

17.663 

(0.025) 

4 

6 

19.671 

(0.039) 

3 

20.581 

(0.023) 

4 

20.416 

(0.039) 

4 

20.215 

(0.035) 

4 

19.591 

(0.062) 

4 

7 


- 

- 

22.328 

(0.102) 

3 

20.957 

(0.014) 

3 

19.984 

(0.018) 

3 

18.869 

(0.046) 

3 

8 

22.416 

(0.099) 

2 

22.312 

(0.098) 

3 

20.987 

(0.031) 

3 

20.103 

(0.009) 

3 

18.957 

(0.028) 

3 

SNF20080522-000 

1 

17.567 

(0.005) 

2 

16.795 

(0.015) 

4 

15.847 

(0.005) 

4 

15.274 

(0.023) 

4 

14.860 

(0.096) 

4 

2 

16.564 

(0.004) 

2 

16.595 

(0.016) 

4 

15.939 

(0.010) 

4 

15.526 

(0.012) 

3 

15.086 

(0.018) 

1 

3 

20.581 

(0.059) 

2 

19.465 

(0.017) 

4 

17.981 

(0.006) 

4 

17.006 

(0.007) 

4 

15.809 

(0.007) 

4 

4 

19.870 

(0.027) 

2 

19.068 

(0.011) 

4 

18.052 

(0.009) 

4 

17.388 

(0.001) 

4 

16.749 

(0.013) 

4 

5 

19.337 

(0.033) 

2 

19.125 

(0.021) 

4 

18.350 

(0.010) 

4 

17.896 

(0.004) 

4 

17.449 

(0.009) 

4 

6 

20.012 

(0.098) 

2 

19.361 

(0.017) 

4 

18.450 

(0.006) 

4 

17.896 

(0.014) 

4 

17.386 

(0.013) 

4 

7 

19.654 

(0.015) 

2 

19.344 

(0.034) 

4 

18.538 

(0.017) 

4 

18.069 

(0.016) 

4 

17.602 

(0.012) 

4 

8 

21.162 

(0.193) 

2 

20.192 

(0.052) 

4 

18.786 

(0.014) 

4 

17.849 

(0.009) 

4 

16.752 

(0.017) 

4 

9 


- 

- 

21.984 

(0.132) 

4 

20.516 

(0.051) 

4 

19.428 

(0.009) 

4 

17.775 

(0.018) 

4 

SNF20080522-001 

1 

15.153 

(0.021) 

1 

15.024 

(0.003) 

2 

14.273 

(0.001) 

2 

13.849 

(0.019) 

1 


- 

- 

2 

15.493 

(0.018) 

2 

15.434 

(0.007) 

3 

14.674 

(0.011) 

4 

14.296 

(0.019) 

1 

13.886 

(0.019) 

1 

3 

16.138 

(0.012) 

2 

15.898 

(0.006) 

3 

15.111 

(0.007) 

4 

14.684 

(0.008) 

4 

14.249 

(0.019) 

1 

4 

16.268 

(0.021) 

2 

16.219 

(0.010) 

3 

15.543 

(0.008) 

4 

15.189 

(0.012) 

4 

14.810 

(0.020) 

3 

5 

18.513 

(0.012) 

2 

17.384 

(0.005) 

3 

16.181 

(0.002) 

4 

15.399 

(0.008) 

4 

14.611 

(0.008) 

2 

6 

17.136 

(0.019) 

2 

16.962 

(0.006) 

3 

16.227 

(0.005) 

4 

15.804 

(0.012) 

4 

15.363 

(0.014) 

4 

7 

17.728 

(0.015) 

2 

17.513 

(0.012) 

3 

16.799 

(0.011) 

4 

16.397 

(0.010) 

4 

16.001 

(0.012) 

4 

8 

18.604 

(0.013) 

2 

17.796 

(0.010) 

3 

16.809 

(0.012) 

4 

16.197 

(0.017) 

4 

15.654 

(0.016) 

4 

9 

17.869 

(0.062) 

2 

18.076 

(0.010) 

3 

17.562 

(0.007) 

4 

17.243 

(0.011) 

4 

16.905 

(0.020) 

4 

10 

19.175 

(0.015) 

2 

18.664 

(0.010) 

3 

17.781 

(0.004) 

4 

17.225 

(0.009) 

4 

16.750 

(0.014) 

4 

11 

18.369 

(0.009) 

2 

18.439 

(0.019) 

3 

17.830 

(0.012) 

4 

17.486 

(0.013) 

4 

17.119 

(0.019) 

4 

12 

20.674 

(0.039) 

2 

19.418 

(0.045) 

3 

18.094 

(0.011) 

4 

17.232 

(0.009) 

4 

16.388 

(0.017) 

4 

13 


- 

- 

20.023 

(0.088) 

3 

18.456 

(0.011) 

4 

17.433 

(0.019) 

4 

16.124 

(0.024) 

4 

14 

19.527 

(0.159) 

2 

19.348 

(0.018) 

3 

18.564 

(0.015) 

4 

18.143 

(0.025) 

4 

17.753 

(0.029) 

4 

15 


- 

- 

19.975 

(0.070) 

3 

18.607 

(0.013) 

4 

17.694 

(0.009) 

4 

16.726 

(0.022) 

4 

16 

18.700 

(0.004) 

2 

19.107 

(0.016) 

3 

18.610 

(0.016) 

4 

18.334 

(0.023) 

4 

18.005 

(0.024) 

4 

17 

18.874 

(0.038) 

2 

19.049 

(0.011) 

3 

18.754 

(0.014) 

4 

18.613 

(0.035) 

4 

18.380 

(0.035) 

4 

18 

21.040 

(0.137) 

2 

20.288 

(0.083) 

3 

18.806 

(0.030) 

4 

17.882 

(0.023) 

4 

16.926 

(0.023) 

4 

19 

19.239 

(0.138) 

2 

19.498 

(0.016) 

3 

19.018 

(0.013) 

4 

18.708 

(0.027) 

4 

18.269 

(0.017) 

3 

20 

21.209 

(0.026) 

2 

20.207 

(0.022) 

3 

19.055 

(0.016) 

4 

18.282 

(0.009) 

4 

17.569 

(0.016) 

4 

21 


- 

- 

21.791 

(0.109) 

3 

19.832 

(0.018) 

4 

18.795 

(0.015) 

4 

17.219 

(0.016) 

4 

SNF20080522-004 

1 

16.881 

(0.001) 

2 

16.539 

(0.006) 

4 

15.707 

(0.013) 

4 

15.216 

(0.010) 

4 

14.790 

(0.086) 

4 

2 

18.114 

(0.020) 

2 

16.996 

(0.008) 

4 

15.867 

(0.014) 

4 

15.181 

(0.013) 

3 

14.607 

(0.019) 

1 
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Table 3. continued. 


No. 

U 

^obs 

B 

^obs 

V 

^obs 

R 

^obs 

/ 

^obs 

3 

17.332 (0.021) 

2 

17.437 (0.006) 

4 

16.804 (0.016) 

4 

16.447 (0.009) 

4 

16.044 (0.007) 

4 

4 

17.313 (0.013) 

2 

17.577 (0.005) 

4 

17.044 (0.020) 

4 

16.708 (0.011) 

4 

16.328 (0.008) 

4 

5 

17.874 (0.052) 

2 

17.911 (0.018) 

4 

17.302 (0.019) 

4 

16.951 (0.003) 

4 

16.579 (0.015) 

4 

6 

18.565 (0.040) 

2 

18.648 (0.026) 

4 

18.066 (0.022) 

4 

17.716 (0.010) 

4 

17.321 (0.004) 

4 

7 

18.768 (0.041) 

2 

18.809 (0.014) 

3 

18.096 (0.013) 

4 

17.703 (0.008) 

4 

17.318 (0.025) 

4 

8 

18.497 (0.054) 

2 

18.758 (0.038) 

4 

18.151 (0.020) 

4 

17.782 (0.006) 

4 

17.354 (0.021) 

4 

9 

19.764 (0.025) 

2 

19.278 (0.033) 

4 

18.444(0.012) 

4 

17.915 (0.017) 

4 

17.405 (0.016) 

4 

10 

18.925 (0.002) 

2 

19.175 (0.014) 

4 

18.668 (0.024) 

4 

18.326 (0.020) 

4 

17.939 (0.019) 

4 

11 

19.639 (0.151) 

2 

19.427 (0.053) 

4 

18.687 (0.020) 

4 

18.219(0.011) 

4 

17.753 (0.026) 

4 

12 

18.793 (0.060) 

2 

19.114(0.032) 

3 

18.694 (0.012) 

4 

18.431 (0.016) 

4 

18.175 (0.037) 

4 

13 

19.789 (0.092) 

2 

19.515 (0.049) 

4 

18.737 (0.023) 

4 

18.297 (0.019) 

4 

17.797 (0.022) 

4 

14 

19.862(0.156) 

2 

19.605 (0.046) 

4 

18.878 (0.026) 

4 

18.424(0.018) 

4 

17.945 (0.006) 

3 

15 

21.046 (0.160) 

2 

20.208 (0.148) 

4 

19.237 (0.043) 

4 

18.578 (0.010) 

4 

17.982 (0.009) 

4 

16 

- 

- 

20.974 (0.077) 

4 

19.544 (0.027) 

4 

18.457 (0.019) 

4 

16.839 (0.005) 

4 

17 

- 

- 

21.699 (0.263) 

3 

20.220 (0.034) 

4 

19.218 (0.010) 

4 

17.857 (0.014) 

4 

SNF20080522-011 

1 

16.565 (0.018) 

1 

16.624 (0.023) 

3 

16.043 (0.017) 

3 

15.667 (0.015) 

3 

15.277 (0.014) 

1 

2 

17.026 (0.019) 

1 

16.825 (0.011) 

3 

16.065 (0.018) 

3 

15.630 (0.017) 

3 

15.194 (0.014) 

1 

3 

18.297 (0.033) 

1 

17.387 (0.006) 

3 

16.358 (0.006) 

3 

15.729 (0.010) 

3 

15.212 (0.014) 

1 

4 

19.371 (0.063) 

1 

18.119(0.112) 

3 

16.840 (0.020) 

3 

15.977 (0.026) 

3 

15.157 (0.014) 

1 

5 

18.045 (0.027) 

1 

17.679 (0.015) 

3 

16.860 (0.006) 

3 

16.371 (0.014) 

3 

15.934 (0.017) 

3 

6 

17.794 (0.025) 

1 

17.841 (0.024) 

3 

17.216 (0.016) 

3 

16.840 (0.014) 

3 

16.465 (0.026) 

3 

7 

17.651 (0.025) 

1 

18.017 (0.025) 

3 

17.544 (0.019) 

3 

17.226 (0.009) 

3 

16.881 (0.028) 

3 

8 

19.910(0.102) 

1 

19.081 (0.047) 

3 

17.782 (0.002) 

3 

16.955 (0.005) 

3 

16.236 (0.012) 

3 

9 

19.422 (0.070) 

1 

18.910(0.011) 

3 

17.996 (0.010) 

3 

17.450 (0.016) 

3 

17.000 (0.025) 

3 

10 

18.491 (0.035) 

1 

18.761 (0.014) 

3 

18.228 (0.018) 

3 

17.889 (0.006) 

3 

17.520 (0.024) 

3 

11 

18.677 (0.041) 

1 

18.932 (0.035) 

3 

18.454 (0.006) 

3 

18.120 (0.008) 

3 

17.747 (0.039) 

3 

12 

20.003 (0.101) 

1 

19.453 (0.012) 

3 

18.535 (0.020) 

3 

17.957 (0.009) 

3 

17.468 (0.030) 

3 

13 

19.470 (0.067) 

1 

19.395 (0.040) 

3 

18.610(0.005) 

3 

18.171 (0.008) 

3 

17.723 (0.020) 

3 

14 

18.925 (0.052) 

1 

19.335 (0.031) 

3 

18.778 (0.009) 

3 

18.432 (0.014) 

3 

18.051 (0.029) 

3 

15 

20.407 (0.161) 

1 

19.925 (0.037) 

3 

18.851 (0.014) 

3 

18.111 (0.008) 

3 

17.486 (0.022) 

3 

16 

20.116(0.135) 

1 

19.907 (0.017) 

3 

19.044(0.031) 

3 

18.503 (0.018) 

3 

18.002 (0.027) 

3 

17 

- 

- 

20.108 (0.005) 

3 

19.057 (0.017) 

3 

18.376 (0.023) 

3 

17.779 (0.024) 

3 

18 

20.185 (0.130) 

1 

19.977 (0.031) 

3 

19.101 (0.007) 

3 

18.590 (0.020) 

3 

18.120 (0.042) 

3 

19 

20.489 (0.168) 

1 

20.081 (0.055) 

3 

19.131 (0.014) 

3 

18.529 (0.016) 

3 

17.924(0.018) 

3 

20 

- 

- 

20.726 (0.065) 

3 

19.356 (0.016) 

3 

18.488 (0.013) 

3 

17.740 (0.030) 

3 

21 

20.057 (0.123) 

1 

20.161 (0.018) 

3 

19.451 (0.020) 

3 

19.016 (0.024) 

3 

18.545 (0.037) 

3 

22 

- 

- 

21.437 (0.167) 

3 

19.971 (0.029) 

3 

18.935 (0.024) 

3 

17.584 (0.033) 

3 

23 

- 

- 

21.387 (0.055) 

3 

20.008 (0.029) 

3 

19.176 (0.018) 

3 

18.328 (0.011) 

3 

SNF20080606-012 

1 

16.230 (0.015) 

1 

16.028 (0.015) 

1 

15.289 (0.012) 

1 

- 

- 

- 

- 

2 

16.481 (0.018) 

2 

16.509 (0.009) 

2 

15.877 (0.010) 

2 

15.521 (0.013) 

1 

- 

- 

3 

17.477 (0.044) 

2 

16.846 (0.013) 

2 

15.882 (0.008) 

2 

15.384 (0.013) 

1 

- 

- 

4 

17.292 (0.015) 

2 

17.233 (0.009) 

2 

16.536 (0.003) 

2 

16.170 (0.002) 

2 

15.733 (0.013) 

1 

5 

17.480 (0.033) 

2 

17.412 (0.005) 

2 

16.748 (0.009) 

2 

16.384 (0.006) 

2 

15.973 (0.013) 

1 

6 

19.395 (0.122) 

2 

18.137 (0.035) 

2 

16.773 (0.009) 

2 

15.920 (0.003) 

2 

- 

- 

7 

19.091 (0.057) 

2 

17.984 (0.014) 

2 

16.832 (0.008) 

2 

16.148 (0.004) 

2 

- 

- 

8 

17.900 (0.012) 

2 

17.822 (0.000) 

2 

17.173 (0.014) 

2 

16.801 (0.011) 

2 

16.382 (0.013) 

1 

9 

19.780 (0.061) 

2 

18.568 (0.001) 

2 

17.304 (0.008) 

2 

16.480 (0.007) 

2 

15.721 (0.013) 

1 

10 

18.034 (0.043) 

2 

17.977 (0.011) 

2 

17.313 (0.011) 

2 

16.951 (0.003) 

2 

16.525 (0.013) 

1 

11 

18.081 (0.022) 

2 

18.049 (0.022) 

2 

17.360 (0.002) 

2 

16.964 (0.010) 

2 

16.514(0.013) 

1 

12 

18.081 (0.022) 

2 

18.049 (0.022) 

2 

17.360 (0.002) 

2 

16.964 (0.010) 

2 

16.514(0.013) 

1 

13 

18.342 (0.043) 

2 

18.226 (0.006) 

2 

17.504 (0.016) 

2 

17.128 (0.004) 

2 

16.706 (0.013) 

1 

14 

19.782 (0.099) 

2 

18.699 (0.011) 

2 

17.654 (0.011) 

2 

17.025 (0.006) 

2 

16.425 (0.013) 

1 

15 

18.405 (0.010) 

2 

18.454 (0.025) 

2 

17.801 (0.009) 

2 

17.422 (0.001) 

2 

16.995 (0.013) 

1 

16 

19.097 (0.084) 

2 

18.716 (0.020) 

2 

17.880 (0.004) 

2 

17.392 (0.004) 

2 

16.882 (0.013) 

1 

17 

18.733 (0.088) 

2 

18.596 (0.009) 

2 

17.887 (0.009) 

2 

17.484 (0.002) 

2 

17.053 (0.013) 

1 

18 

18.779 (0.054) 

2 

18.739 (0.008) 

2 

18.072 (0.014) 

2 

17.707 (0.012) 

2 

17.290 (0.013) 

1 

19 

18.804(0.070) 

2 

18.842 (0.029) 

2 

18.201 (0.005) 

2 

17.834 (0.001) 

2 

17.426 (0.013) 

1 

20 

18.802(0.078) 

2 

18.896 (0.013) 

2 

18.261 (0.006) 

2 

17.913 (0.008) 

2 

17.475 (0.013) 

1 
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Table 3. continued. 


No. 

U 

^obs 

B 

^obs 

V 

^obs 

R 

^obs 

/ 

^obs 

21 

19.810(0.161) 

2 

19.290 (0.052) 

2 

18.377 (0.012) 

2 

17.849 (0.004) 

2 

17.345 (0.013) 

1 

22 

18.918 (0.083) 

2 

18.998 (0.001) 

2 

18.386 (0.015) 

2 

18.028 (0.008) 

2 

17.625 (0.013) 

1 

23 

18.902 (0.001) 

2 

19.146 (0.015) 

2 

18.663 (0.002) 

2 

18.367 (0.000) 

2 

17.977 (0.014) 

1 

24 

19.064 (0.137) 

2 

19.246 (0.017) 

2 

18.718 (0.021) 

2 

18.399 (0.009) 

2 

17.989 (0.014) 

1 

25 

18.947 (0.037) 

2 

19.233 (0.026) 

2 

18.722 (0.000) 

2 

18.402 (0.003) 

2 

17.974 (0.014) 

1 

26 

20.340 (0.195) 

2 

19.769 (0.016) 

2 

18.838 (0.015) 

2 

18.299 (0.000) 

2 

17.721 (0.014) 

1 

27 

21.102 (0.027) 

2 

20.042 (0.007) 

2 

18.883 (0.020) 

2 

18.103 (0.012) 

2 

17.387 (0.014) 

1 

28 

- 

- 

19.899 (0.021) 

2 

19.044 (0.033) 

2 

18.533 (0.001) 

2 

18.015 (0.014) 

1 

29 

20.117(0.058) 

2 

19.866 (0.003) 

2 

19.047 (0.003) 

2 

18.564 (0.007) 

2 

18.075 (0.014) 

1 

30 

- 

- 

20.015 (0.011) 

2 

19.097 (0.024) 

2 

18.611 (0.003) 

2 

18.079 (0.014) 

1 

31 

19.506 (0.119) 

2 

19.712(0.022) 

2 

19.160 (0.009) 

2 

18.808 (0.010) 

2 

18.411 (0.015) 

1 

32 

19.473 (0.028) 

2 

19.727 (0.013) 

2 

19.174 (0.021) 

2 

18.841 (0.020) 

2 

18.431 (0.015) 

1 

33 

- 

- 

20.668 (0.035) 

2 

19.251 (0.007) 

2 

18.349 (0.004) 

2 

17.393 (0.014) 

1 

34 

19.758 (0.012) 

2 

19.967 (0.052) 

2 

19.279 (0.018) 

2 

18.823 (0.008) 

2 

18.325 (0.015) 

1 

35 

- 

- 

20.053 (0.026) 

2 

19.355 (0.024) 

2 

18.895 (0.006) 

2 

18.384 (0.015) 

1 

36 

19.604 (0.063) 

2 

19.936 (0.026) 

2 

19.388 (0.000) 

2 

19.056 (0.012) 

2 

18.691 (0.016) 

1 

37 

- 

- 

20.983 (0.037) 

2 

19.440 (0.009) 

2 

18.375 (0.001) 

2 

16.661 (0.015) 

1 

38 

20.071 (0.160) 

2 

20.125 (0.007) 

2 

19.545 (0.014) 

2 

19.255 (0.032) 

2 

18.843 (0.017) 

1 

39 

- 

- 

21.407 (0.019) 

2 

19.666 (0.007) 

2 

18.516(0.008) 

2 

16.701 (0.015) 

1 

40 

- 

- 

21.052 (0.047) 

2 

19.709 (0.029) 

2 

18.901 (0.020) 

2 

18.061 (0.015) 

1 

41 

20.397 (0.099) 

2 

20.544(0.081) 

2 

19.769 (0.003) 

2 

19.332 (0.022) 

2 

18.865 (0.017) 

1 

42 

19.926 (0.079) 

2 

20.224 (0.018) 

2 

19.772 (0.010) 

2 

19.448 (0.013) 

2 

19.055 (0.018) 

1 

43 

- 

- 

21.079 (0.022) 

2 

19.799 (0.033) 

2 

18.992 (0.015) 

2 

18.224 (0.015) 

1 

44 

- 

- 

21.129 (0.081) 

2 

19.903 (0.006) 

2 

19.119(0.001) 

2 

18.393 (0.016) 

1 

45 

- 

- 

21.561 (0.102) 

2 

19.992 (0.014) 

2 

19.008 (0.026) 

2 

17.680 (0.015) 

1 

46 

- 

- 

21.395 (0.056) 

2 

20.017 (0.008) 

2 

19.188 (0.000) 

2 

18.335 (0.016) 

1 

47 

- 

- 

20.554 (0.012) 

2 

20.039 (0.038) 

2 

19.746 (0.018) 

2 

19.324 (0.020) 

1 

48 

- 

- 

21.196 (0.027) 

2 

20.238 (0.005) 

2 

19.589 (0.004) 

2 

18.948 (0.018) 

1 

49 

- 

- 

21.767 (0.011) 

2 

20.489 (0.015) 

2 

19.617 (0.025) 

2 

18.769 (0.017) 

1 

50 

22.084 (0.200) 

2 

22.137 (0.034) 

2 

20.587 (0.010) 

2 

19.597 (0.005) 

2 

18.128 (0.016) 

1 

51 

- 

- 

22.048 (0.051) 

2 

20.619 (0.004) 

2 

19.661 (0.015) 

2 

18.401 (0.016) 

1 

52 

- 

- 

21.813 (0.243) 

2 

20.692 (0.002) 

2 

19.993 (0.018) 

2 

19.276 (0.021) 

1 

SNF20080610-003 

1 

15.941 (0.014) 

1 

15.862 (0.014) 

1 

15.215 (0.012) 

1 

14.846 (0.012) 

1 

14.582 (0.012) 

1 

2 

17.660 (0.022) 

1 

16.605 (0.017) 

1 

15.506 (0.012) 

1 

14.842 (0.013) 

1 

- 

- 

3 

16.781 (0.016) 

1 

16.330 (0.015) 

1 

15.551 (0.012) 

1 

15.112(0.012) 

1 

- 

- 

4 

16.343 (0.014) 

1 

16.225 (0.014) 

1 

15.552 (0.012) 

1 

15.169 (0.012) 

1 

14.753 (0.012) 

1 

5 

18.616 (0.030) 

2 

17.283 (0.018) 

2 

15.975 (0.007) 

2 

15.157 (0.013) 

1 

- 

- 

6 

16.493 (0.037) 

2 

16.560 (0.016) 

2 

16.002 (0.001) 

2 

15.685 (0.012) 

1 

15.327 (0.012) 

1 

7 

16.549 (0.015) 

1 

16.698 (0.013) 

1 

16.161 (0.012) 

1 

15.817(0.013) 

1 

15.438 (0.012) 

1 

8 

16.910(0.020) 

2 

16.946 (0.011) 

2 

16.324 (0.008) 

2 

15.973 (0.002) 

2 

15.608 (0.045) 

2 

9 

17.444 (0.040) 

2 

17.210 (0.027) 

2 

16.460 (0.000) 

2 

16.065 (0.002) 

2 

15.643 (0.012) 

1 

10 

16.742 (0.031) 

2 

17.022 (0.031) 

2 

16.549 (0.003) 

2 

16.261 (0.000) 

2 

15.906 (0.013) 

2 

11 

17.815 (0.030) 

2 

17.476 (0.013) 

2 

16.704 (0.007) 

2 

16.274 (0.007) 

2 

15.870 (0.012) 

1 

12 

17.559 (0.019) 

2 

17.474 (0.020) 

2 

16.781 (0.006) 

2 

16.383 (0.013) 

2 

15.939 (0.006) 

2 

13 

17.987 (0.024) 

2 

18.151 (0.006) 

2 

17.495 (0.004) 

2 

17.086 (0.013) 

2 

16.610(0.012) 

2 

14 

17.992 (0.007) 

2 

18.166 (0.015) 

2 

17.632 (0.003) 

2 

17.318 (0.011) 

2 

16.953 (0.012) 

2 

15 

18.471 (0.030) 

2 

18.349 (0.020) 

2 

17.658 (0.008) 

2 

17.278 (0.005) 

2 

16.872 (0.005) 

2 

16 

20.765 (0.141) 

2 

19.453 (0.011) 

2 

18.046 (0.000) 

2 

17.170 (0.009) 

2 

16.267 (0.010) 

2 

17 

20.008 (0.085) 

2 

19.169 (0.038) 

2 

18.173 (0.011) 

2 

17.599 (0.001) 

2 

17.090 (0.016) 

2 

18 

19.297 (0.065) 

2 

19.033 (0.008) 

2 

18.275 (0.011) 

2 

17.848 (0.003) 

2 

17.423 (0.013) 

2 

19 

20.263 (0.010) 

2 

19.369 (0.070) 

2 

18.322 (0.010) 

2 

17.637 (0.001) 

2 

17.036 (0.013) 

2 

20 

19.827 (0.043) 

2 

19.348 (0.013) 

2 

18.493 (0.001) 

2 

18.019 (0.014) 

2 

17.553 (0.007) 

2 

21 

- 

- 

20.202 (0.075) 

2 

18.665 (0.020) 

2 

17.697 (0.010) 

2 

16.532 (0.023) 

2 

22 

18.960 (0.005) 

2 

19.174 (0.056) 

2 

18.704(0.003) 

2 

18.406 (0.016) 

2 

18.054 (0.025) 

2 

23 

19.255 (0.016) 

2 

19.504 (0.016) 

2 

19.059 (0.006) 

2 

18.781 (0.001) 

2 

18.456 (0.008) 

2 

24 

20.067 (0.051) 

2 

19.904 (0.017) 

2 

19.144(0.006) 

2 

18.757 (0.029) 

2 

18.283 (0.002) 

2 

25 

- 

- 

20.663 (0.012) 

2 

19.160 (0.012) 

2 

18.187 (0.003) 

2 

16.992 (0.029) 

2 

26 

21.732 (0.032) 

2 

20.695 (0.022) 

2 

19.174 (0.007) 

2 

18.188 (0.001) 

2 

16.892 (0.002) 

2 

27 

- 

- 

20.227 (0.004) 

2 

19.347 (0.007) 

2 

18.816(0.009) 

2 

18.302 (0.008) 

2 
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Table 3. continued. 


No. 

U 

^obs 

B 

^ obs 

V 

^ obs 

R 

^ obs 

/ 

^obs 

28 

- 

- 

20.608 (0.006) 

2 

19.347 (0.026) 

2 

18.548 (0.007) 

2 

17.791 (0.007) 

2 

29 

19.510(0.004) 

2 

19.815 (0.041) 

2 

19.349 (0.011) 

2 

19.026 (0.011) 

2 

18.682 (0.001) 

2 

30 

- 

- 

21.039 (0.027) 

2 

19.386 (0.006) 

2 

18.351 (0.001) 

2 

16.825 (0.012) 

2 

31 

18.971 (0.034) 

2 

19.954 (0.020) 

2 

19.552 (0.010) 

2 

19.399 (0.007) 

2 

19.195 (0.017) 

2 

32 

21.493 (0.188) 

2 

20.610(0.043) 

2 

19.595 (0.013) 

2 

18.906 (0.006) 

2 

18.321 (0.027) 

2 

33 

21.461 (0.194) 

2 

20.654 (0.069) 

2 

19.609 (0.007) 

2 

18.961 (0.013) 

2 

18.301 (0.009) 

2 

34 

- 

- 

21.073 (0.057) 

2 

19.644 (0.014) 

2 

18.707 (0.002) 

2 

17.530 (0.026) 

2 

35 

20.202 (0.148) 

2 

20.334 (0.048) 

2 

19.658 (0.001) 

2 

19.268 (0.046) 

2 

18.798 (0.015) 

2 

36 

- 

- 

20.694 (0.005) 

2 

19.667 (0.015) 

2 

19.036 (0.003) 

2 

18.472 (0.032) 

2 

37 

- 

- 

21.194 (0.052) 

2 

19.720 (0.021) 

2 

18.705 (0.014) 

2 

17.205 (0.005) 

2 

38 

- 

- 

21.016(0.000) 

2 

19.906 (0.001) 

2 

19.162 (0.017) 

2 

18.502 (0.013) 

2 

39 

- 

- 

21.765 (0.076) 

2 

20.178 (0.030) 

2 

19.167 (0.008) 

2 

17.691 (0.018) 

2 

40 

20.458 (0.038) 

2 

20.739 (0.054) 

2 

20.264 (0.009) 

2 

19.973 (0.006) 

2 

19.649 (0.006) 

2 

41 

- 

- 

22.000 (0.049) 

2 

20.326 (0.008) 

2 

19.349 (0.003) 

2 

17.902 (0.011) 

2 

42 

- 

- 

21.666 (0.021) 

2 

20.369 (0.007) 

2 

19.560 (0.039) 

2 

18.833 (0.003) 

2 

43 

- 

- 

22.120 (0.004) 

2 

20.610 (0.008) 

2 

19.695 (0.023) 

2 

18.563 (0.008) 

2 



